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ABSTRACT 


Displacement theories are developed for a variety of structures with the goal of providing real- 
time shape predictions for aerospace vehicles during flight. These theories are initially developed 
for a cantilever beam to predict the deformed shapes of the Helios flying wing. The main structural 
configuration of the Helios wing is a cantilever wing tubular spar subjected to bending, torsion, 
and combined bending and torsion loading. The displacement equations that are formulated are 
expressed in terms of strains measured at multiple sensing stations equally spaced on the surface 
of the wing spar. The bending and torsion strain data are to be input to the displacement equations 
for calculations of slopes, deflections, and cross-sectional twist angles of the wing spar at the 
strain-sensing stations. Displacement theories for other structures, such as tapered cantilever 
beams, two-point supported beams, wing boxes, and plates also are developed. The accuracy of 
the displacement theories is successfully validated by finite-element analysis and classical beam 
theory using input-strains generated by finite-element analysis. The displacement equations and 
associated strain-sensing system (such as fiber optic sensors) create a powerful means for in-flight 
deformation monitoring of aerospace structures. Because local strains are used as inputs for the 
structural deformation predictions described in this report, this method serves multiple purposes 
for structural shape sensing, loads monitoring, and structural health monitoring. Ultimately, the 
calculated displacement data can be visually displayed to the ground-based pilot or used as input 
to the control system to actively control the shape of structures during flight. The displacement 
theories could also be easily applied to calculate the deformed shapes of a variety of high-aspect- 
ratio flexible structures, aircraft wing rotorcraft vehicles, space vehicles, wind turbine blades (for 
the alternative energy community), space-based structures (such as booms and antennae), and 
long-span civil structures (such as bridges and dams). 

NOMENCLATURE 

A cross-sectional area of uniform cantilever beam, in 2 

A1 aluminum 

a mean radius of 4-ply composite tube wall and aluminum tube wall, or mean radius 

of Helios 31 -ply reinforced spar cap wall, in. 

c outer radius of tubular beam, or half depth of beam, in. 

c) distance from the neutral axis to the surface strain-sensing station, x t , in. 

c n value of c, at beam tip (i = n ), in. 

c 0 value of c, at beam root (i = 0), in. 

d chord-wise distance between front and rear sensor lines for cantilever wing 

box, in. 

deg degree 

E Young modulus, lb/in 2 

E l Young modulus of lamina in fiber direction, lb/in 2 

E T Young modulus of lamina in direction transverse to fiber direction, lb/in 2 



finite-element span- wise length, in. 
shear modulus, lb/in 2 
shear modulus of lamina, lb/in 2 
thickness of plate, in. 

moment of inertia of tube cross section, in 4 

= 0, 1,2, 3, n, strain-sensing station identification number 

polar moment of inertia, in 4 
index 

length of cantilever tube, in. 
moment, in-lb 

moment at strain-sensing station, x t , in-lb 

number associated with the last span- wise strain- sensing station 

applied load, lb 

Structural Performance And Resizing 
twisting moment, in-lb 

twisting moment at the distortion-sensing station, x t , in-lb 
thickness of 4-ply composite wall (or aluminum wall), in. 
thickness of 31-ply composite reinforced wall region, in. 
local shear force at strain-sensing station, v, , lb 
width of wing box, in. 

Cartesian coordinates, in. 

deflection in y-dircction at any point, a;, between two strain-sensing stations, in. 

i - th strain-sensing station, or its axial coordinate, x= x t , in. 

deflection in y-dircction at strain-sensing station, x t , in. 

deflection of two-point supported beam at strain- sensing station, x t , in. 

shear-induced deflection, in. 

deflection in y-dircction of the wing box rear sensing line at strain-sensing station, 

x i , in. 

surface distortion angle (shear strain) at distortion strain-sensing station, x t , rad 
deflection shear strain at strain- sensing station, x i , rad 

= l/n, distance between two adjacent strain-sensing stations (equally spaced), in. 
displacements at two ends of finite element, in. 
bending strain, in/in. 

bending strain at strain-sensing station, x t , due to bending only, in/in. 

bending strain of the wing box rear sensing line at strain-sensing station, , due to 
bending only, in/in. 



true bending strain at strain-sensing station, x t , due to bending and torsion, in/in. 

principal tensile strain at distortion-sensing station, x j , in the 45-deg helical 
direction, in/in. 

beam slope at any point, x, between two strain-sensing stations, x i-\ - x - •*'/ , rad 
or deg 

beam slope at strain-sensing station, x j , rad or deg 
Poisson’s ratio of lamina 
Poisson’s ratio 

bending stress of outermost fiber, lb/in 2 

principal tensile stress at distortion-sensing station, x t , in the 45-deg helical 
direction, lb/in 2 

maximum shear stress at distortion-sensing station, v, , lb/in 2 
cross-sectional twist angle at strain-sensing station, x t , rad or deg 
outboard side of strain-sensing station, x i 
inboard side of strain-sensing station, x t 

INTRODUCTION 

The ultralightweight Helios flying wing (fig. 1; see Appendix A for details) has a wingspan 
of 247 ft. The Helios wingspan is longer than the wingspan of the U.S. Air Force C-5 military 
transport (Lockheed Martin Aeronautics Company, Marietta, Georgia), which is 222 ft, and that of 
the Boeing 747 commercial jetliner (Boeing Company, Seattle, Washington), which ranges from 
195 to 215 ft depending on the model. The major load-carrying wing main spar is a carbon fiber 
composite tube (radius of 4 in.), which is thicker on the top and bottom (spar caps) to improve 
bending stiffness. The spar is wrapped with Nomex® and Kevlar® (both registered trademarks of 
E.I. duPont de Nemours and Company, Wilmington, Delaware) for additional strength. 

Because of an excessively long wingspan and high flexibility, the Helios wing regularly 
experienced large deformations during flight, with the wingtip deflections reaching as high as 
40 ft. This kind of deformed shape is similar to a 72-deg circular arc with a 197-ft radius. 

The Helios wing recently experienced a midair breakup while flying at an altitude of 
approximately 3000 ft. The primary cause of the mishap is believed to have resulted from undamped 
pitch oscillations of the highly deformed wing. 

To avoid future mishaps, a method must be developed in which the deformations of highly 
flexible structures can be visually monitored during flight. This way, if the wingtip deflections 
approach the design limit, the ground-based pilot could initiate an emergency maneuver to avoid 
further deformations. 
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This report discusses the formulation of displacement theories for uniform and tapered 
cantilever beams under bending, torsion, and combined bending and torsion. The formulated 
displacement equations are expressed in terms of strains at multiple strain sensors attached to the 
wing spar surface and can be used to calculate wing deflections and twist angles during flight. Next, 
the displacement theories for other structural configurations, such as two-point supported beams, 
wing boxes, and plates also are developed. Finally, the accuracy of the displacement theories is 
validated with finite-element analysis and classical beam theory using input-strains generated by 
finite-element analysis. 

The displacement equations and onboard strain-sensing system could be a powerful tool for 
in-flight monitoring of wing deformations by the ground-based pilot for maintaining safe flights 
of the unmanned flying wings. In addition, the real-time wing shape could be input to the aircraft 
control system for aeroelastic wing shape control. 

IN-FLIGHT DEFORMATION MONITORING 

In-flight wing deformation monitoring of the Helios ultralightweight flying wing is particularly 
difficult because of weight restrictions and highly flexible nature of the aircraft. Several methods 
for in-flight monitoring of wing deflections exist for use on conventional aircraft. 

One method is the electro-optical flight deflection measurement system, which is composed of 
onboard cameras and several wing-mounted targets (ref. 1). This system provides wing displacement 
information during the flight but is too heavy for lightweight flying wing applications. 

Another in-flight deflection measurement method is the use of conventional strain gages to 
collect local strains for wing deflection information. Numerous strain-sensing stations are required 
to capture the higher displacement modes of the flexible wings. When multiple strain-sensing 
stations are used, the weight of the strain gage lead wires alone is too heavy and impractical for 
most weight conscious lightweight flying wing aircraft. 

Fiber optic sensors offer the most attractive alternative to conventional strain gages. They are 
lightweight, fine, and flexible filaments (approximately the size of human hairs) and can be highly 
multiplexed at desired sensing intervals. 

THEORETICAL APPROACH 

If multiple strain sensors are installed at discrete sensing stations on a cantilever wing spar, 
strain sensor data can be used to calculate the deflections and twists of the wing spar during flight. 
Classical beam theory can be used to develop theoretical slope, deflection, and cross-sectional 
twist angle equations for uniform cantilever beams. These displacement equations are written in 
terms of strains for the calculations of the deformed shape of the beam. In the present report, the 
input strains are calculated from the displacement outputs of finite-element stress analysis. The 
accuracy of the newly developed displacement equations (for bending and torsion) is then verified 
by classical cantilever beam theory. 
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With slight modifications, the displacement equations developed for the uniform cantilever 
beam could be applied to predict the deformed shapes of two-point supported beams (simple 
beams) and plates. The displacement equations also can be modified to formulate the displacement 
equations for tapered, slightly tapered, and stepwise tapered cantilever beams, and for tapered 
cantilever wing boxes. 


UNIFORM CANTILEVER BEAMS 

Basic displacement equations are developed for a uniform cantilever beam (for example, 
the Helios wing tubular spar) subjected to bending, torsion, and combined bending and torsion. 
Although the highly flexible Helios wing can undergo large deformation (wingtip deflections of 
as much as 40 ft), the local strains on the wing spar surface remain small. Therefore, the classical 
small strain theory can be used to derive the slope and displacement equations for the uniform 
cantilever beam. 


Bending 

In this section, displacement equations for bending only are formulated. These equations 
describe the local slopes and deflections at points along the strain- sensing line on the bottom 
surface of a uniform cantilever beam of arbitrary cross-sectional shape (for example, a tubular 
Helios spar). 

Moment- Strain Relationship 

The classical beam differential equation (elastic curvature of the uniform beam) is given by 
(refs. 2, 3) 


d 2 y _ M (x) 
dx 2 El 

in which y is the vertical displacement, x is the span- wise coordinate, Mix) is the bending moment, 
E is the Young modulus, and 1 is the moment of inertia. The bending moment, M(x), in equation (1) 
may be related to the associated bending strain, e(x ) , at the bottom or top fiber (for example, a Helios 
tubular spar generatrix, which is a straight line on the tubular surface parallel to the longitudinal 
axis, x, of the tube) through the outermost fiber bending stress equation, <j(x) = \M(x)c]/l 
(refs. 2, 3): 


e(x) = 


o(x) 

E 


M(x)c 

El 


o(x) = 


M(x)c 

I 


( 2 ) 


In equation (2), c is the beam half depth (for example, the outer radius of the Helios tubular spar). 
In light of equation (2), beam bending equation (1) can be written in terms of e(x) as 
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( 3 ) 


d y _ £(x) 
dx 2 c 

Note that under the strain formulation, beam bending equation (3) contains only the beam 
geometrical parameter, half depth c, and that the flexural rigidity term El is eliminated. 

Consider the case of a highly flexible, extremely long Helios wing tubular spar that carries 
concentrated weight (solar powered motors, pods, hydrogen storage tanks, and so forth) at different 
wing stations. Because the wing is subjected to aerodynamic lift forces during flight, the wing spar 
bending moment, M(x), could be a complicated nonlinear function of the span- wise coordinate, x. 

Consider the uniform cantilever tubular spar with an outer radius, c, and length, l, as shown in 

figure 2. Let the surface bending strain sensors be installed at n + 1 (0, 1, 2, n) number of equally 

spaced strain- sensing stations along the bottom or top generatrix of the spar (tig. 2). Although the 
bending strain is zero at the spar tip, the nth bending strain sensor is installed there because of 
mathematical convenience in the derivation of the spar tip slope and deflection equations. 

The bending moment, M(x), can be assumed to be a piecewise-linear function along the 
axial coordinate, x. Namely, in the region, x, , < x < x,-, between any two adjacent strain-sensing 
stations, { x,-_, x ; }, M(x) is considered as a linear function of (x - x, , ) as 


M(x) = M,.., - (M._, - M,) X *'~ 1 ; x i _ l <x<x i (4) 

A l 

in which {M ( . are the bending moments at the two adjacent strain-sensing stations, (x ; , , x ( } , 
respectively, and A / ( = x t - x ( _, = l/n ) is the axial distance between the two adjacent strain-sensing 
stations, {x ( . p x ( } . In light of equation (2), the local moment equation (4) can be written in terms 
of the local bending strain, £(x) , for the region, x ( , < x < x,. , as 


£(x) = -e,) X Xl ~' ; x ; _j < x < x ; ; A / = - (5) 

A / n 

in which {e i , £ ( } are the bending strains measured at the two adjacent strain- sensing stations, 
{ x,._j , x ( . }, respectively. 

Slope Equations 

The slope, tan0(x) , of the uniform beam in the region, x ( _, < x < x ; , between the two adjacent 
strain-sensing stations, {x ( _,,x, j , can be obtained by integrating equation (3), with the constant of 
integration determined by enforcing the continuity of slope at the inboard adjacent strain-sensing 
station, x,_, , as 
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(6) 


tan(?(x) = [ 

J x i 


d 2 y 
dx 2 


Slope increment 


dx + tan 61, j 

Slope at x,_j 


£(X) 


dx + tan 0, 


i - 1 


A- 1 < X < X,. 


in which tan0 ( _, (constant of integration) is the slope at the inboard adjacent strain-sensing station, 
x,_| . By substituting strain equation (5) into slope equation (6) and carrying out the integration, 
one obtains the slope, tan 61 = [ tan 0( x, ) ] , at the strain-sensing station, x ; , as (see Appendix B for 
details of integration) (ref. 4) 

A l , , 

tan 61. = —(£,_, + £,.) + tan0,._, ; (i = 1, 2, 3, n) (7) 


When i'=l, the slope, tan0 [ , = tan 0 O , will be the slope at the left fixed end, which is zero 
(tand 0 = 0) for the cantilever beam. If the reference left end is simply supported (as in the two-point 
supported beam case), then the slope, tan 0 t] , is nonzero ( tan 0 <t ^ 0) . 

Writing out equation (7) for different strain-sensing stations, (i= 1, 2, 3, , n), results in 


A/ 

— I 

2c 

A/ 

2c* 

A/ 

: 1 

2c 


A / 
2c 
A l_ 
2c 


e 0 + + £ n 


7=1 


; ( tan 0 <t = 0 at built-in end) 

(8-a) 

M 

— (£ 0 + 2f, +£ 2 ) 
2c 

(8-b) 

A/r 1 

— [£ 0 + 2(£j +£ 2 ) + £ 3 J 

(8-c) 

Mr , 

: ~[ £ 0 + 2 (£| +e 2 + £ 3 ) + e 4 J 

2 c 

(8-d) 

A/r 1 

-i - ~l e 0 + 2 (£| +e 2 +£3 + ....+ £„_i)+£„J 
2c 


(8-e) 


The primary objective is to calculate tan 61 (i= 1, 2, 3 , n) for all strain-sensing stations. Because 

the last station, x n , is located at the beam tip, the slope, tan 61, , will give the beam tip slope. 
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Deflection Equations 

The deflection, y(x) , of the uniform beam in the region, x j _ l < x< x t , between the two 
adjacent strain-sensing stations, , can be obtained by integrating slope equation (6) with 

the constant of integration determined by enforcing the continuity of deflection at the inboard 
adjacent strain-sensing station, x ( _, , as 


y(x) = f tan 0(x)dx + 

J X; , 


Integration of slope 


Deflection 

at*.-. 



Deflection increment 


tan0 ; , dx + y t , 

v ' Deflection 

Deflection at x at 

due to tan0,_j 




(9) 


in which y { _ l (integration constant) is the deflection at the inboard adjacent strain- sensing station, 
x t _ t . Carrying out the integration of equation (9) in light of strain equation (5), the deflection, 
y. = [y(x ( )J , at the strain-sensing station, x t , is obtained as (see Appendix B for details of 
integration) (ref. 4) 


y< 


(A/) 2 

6c 


(2 £,_! + £, ) + y,_! + M tan0,._ l ; (i = 1, 2, 3, n) 


( 10 ) 


When i = 1, y 0 = tan 0 O = 0 at the fixed end. Writing out deflection equation (10) for different 
strain-sensing stations (i=l, 2, 3, ...., n ), one obtains 

(A/) 2 , x (A if . x 

y 1 = ( 2e 0 + ej + y 0 + A/tan0 o = (2e 0 +e 1 ) ; y 0 = tan 0 O = 0 at built-in end (10-a) 

6c ' 6c 


(A If , x (A If . N 

_y 2 — — y — (2£| + £ 2 ) + V; + A/ tan0j — (2£ 0 T 3£j + £ 2 j + A l tan0j 


6c 


6c 


(10-b) 


y 3 — — — — (2£ 2 + £ 3 ^ + y 2 + A l tan 0 2 — — — — j^2£ 0 + 3 ^£j + £^^ + £^J + A l ( tan 0 j + tan 0 7 ^ (10-c) 


y 4 


(A If 

6c 


(2£ 3 + £ 4 ) + y 3 + A l tan0, — ^2£ 0 T 3^£ 3 + £, + £ 3 ) + £ 4 J 


6c 


(10-d) 


+ A / ( tan 0, + tan 0 2 + tan 0, ) 


y n — 1” 2£ 0 + 3 f £j + £ 2 + £ 3 + + £„_[ ) + £„! + A / (tan0, + tan0 2 + tan0 3 + + tan0 (! _j 1 

6c 1 ~ ' 

(10-e) 


(MY 


6c 


n - 1 


2e 0 + + e n 


n - 1 

+ A/^T tan 0, 

i=i 
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In light of slope equation (8), the summation of tan (9 of the last term in equation (10-e) can 
be written out in terms of strains, £ . , as 


»- | A / 

Xtan0, = — {(«- lR + [2(«- 1)- ljf, + [2(n- 2)- l]e 2 +\2 (n- 3)- l]e, + ....+ e n _j} 
2c 


i— 1 


A/ 

2c 


("-iR+E[ 2 ( n_ 0-iR 


( 11 ) 


1=1 


In light of equation (11), bending deflection equation (10) takes on the following form, expressed 
exclusively in terms of the bending strains, £ : 


(A If , N 

v, = (2e 0 + gj J ; (y 0 = tan 0 O = 0 at the built-in end) 

6c 


(12-a) 


(A/) 2 , x 

y 2 — (5e 0 + 6e t +e 2 ) 
6c 


(12-b) 


^3 


(A/) 2 

6c 


[8e 0 + 6(2e, +£ 2 ) + e 3 ] 


(12-e) 


(A/) 2 r / 3-1 

y 4 [ll£ 0 +6(3£, +2£ 2 +£ 3 ) + £ 4 J 


(12-d) 


(Air 

6c 


■|(3n — 1)£ 0 + 6^(/r — 1)£] + (n ~ 2)£ 9 + (n — 3)£ 3 + + \n — (n — l)]£ w _j + £ n j - 


(A/) 2 

6c 


n — 1 

(3 n - 1)£ 0 + 6 J, (n - i)£ f + £„ 
/=1 


(12-e) 


The deformed shape of the tubular wing spar can be constructed by plotting each deflection, 
y. ( i = 1, 2, 3...., n), at each strain- sensing station, x ( . The shape could then be displayed on the 
screen so that the ground-based pilot can view the real-time wing deformation status. 

The contribution of deflection caused by transverse shear effect usually is very small for a 
long cantilever beam. If the transverse shear effect is desired, however, an additional transverse 
shear effect term must be added to equation (12) as shown in Appendix C. 
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Torsion 


In this section, displacement equations are developed for a cantilever tubular beam under 
torsion. The local beam cross-sectional twist angle is written in terms of surface distortion angle at 
the inboard distortion-sensing stations. 

Torque-Twist Relationship 

Consider a cantilever tubular spar, with a length of l and n + 1 (0, 1, 2, ... . n) number of 
equally spaced distortion strain-sensing stations (fig. 3), subjected to varying span-wise twisting 
moments in addition to bending (fig. 2). Let 7j be the twisting moment at the distortion-sensing 
station, x i , and let y ; be the associated local surface distortion angle (shear strain) at station x i 
(fig. 3). The torque-distortion relationship can then be written as (refs.l, 2) 



(13) 


in which G is the shear modulus, J is the polar moment of inertia, and c is the outer radius of the 
wing spar. The surface distortion angle, y i (shear strain), in equation (13) is obtained from the 
distortion sensors as 


Yi = ' 


(t ). 

v max s i 


= —£■ = 2(1 + v)£- 


(14) 


in which (r max ), is the maximum shear stress at the distortion-sensing station, x t , and of 
(=(T max ),) and £■' are the principal tensile stress and strain, respectively, in the 45-deg helical 
direction at the distortion-sensing station, v, (fig. 3). 


Cross-Sectional Twists 

In light of figure 3, the total cross-sectional twist angle, 0 ; , at the strain- sensing station, x i , 

can be obtained by summing up the partial cross-sectional twist angles, — y . , at all the inboard 
sensing stations, j (= 1, 2, 3 , /— 1 ) , as (fig. 3) c 


0i = 


A l 

— 7o 
c 


(15-a) 


02 


— (7 0 + 7i) 
c 


(15-b) 


03 


— (7 0 + 7i+7 2 ) 

c 


(15-c) 
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(15-d) 



t = — (Yo + r i + Yi + r 3 + r 4 + •••• + r„_ i) ( 15 -e) 

c 

Notice that in the expression of the tip cross-sectional twist angle, </) n (eq. (15-e)), the tip distortion 
term, y n , does not appear. Therefore, installation of the tip distortion sensor for y n is not required 
(fig. 3). 

Equation (15) can be applied to a rectangular cantilever beam (for example, a wing box) if 
c represents the half width of the wing box. For the wing box, the surface distortion sensors can 
be installed on the vertical or horizontal surface. A simpler way to obtain the cross-sectional twist 
angle for the wing box, however, is to use two parallel bending strain-sensing lines without the 
need to install the distortion sensors (see the section entitled “Other Applications: Wing Boxes”). 

Combined Bending and Torsion 

In this section, displacement equations are developed for a tubular cantilever spar under 
combined bending and torsion. The calculations of true bending strains are presented, from which 
modified slope and modified deflection equations are established. 

True Bending Strain 

As shown in figure 3, if the tubular spar is subjected to both bending and torsion, the bending 
strain-sensing station, x t (for measuring the bending strain, e ; ), will move tangentially by a 
cross-sectional twist angle, 0, , and the bending strain-sensing axis will be tilted by a local surface 
distortion angle, y f . Thus, the strain output, e ; , obtained at the strain-sensing station, x t , no longer 
gives the true bending strain. To obtain the true bending strain, e j , the value of the measured strain, 
e j , must be corrected by using the bending strain correction equation 



cos 0, cos y ( . 


(16) 


in which the surface distortion angle, y, , and the cross-sectional twist angle, (f ) ! , are calculated from 
equations (14) and (15), respectively. With the effect of torsion, slope equation (8) and deflection 
equation (12) for the bending must be modified by replacing e, with e, . 
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Modified Slope Equation 


By replacing the bending strain, £ i , in slope equation (8) with the true bending strain, £ ( . , one 
obtains the following modified slope equation: 


A/ _ _ 

tan0, = —(£ 0 + e i) 
2c 


(17-a) 


A/ _ _ A/ _ _ _ 

tan Oj — — (£| + £^ ) + tan 0 j — — ( £ 0 + 2e 1 + £ 0 ) 
2c 2c 


(17-b) 


A / _ _ A/ r _ _ _ _ n 

tan0, (e 2 + £ 3 ) + tan0 2 [e 0 + 2(e, + £ 2 ) + £ 3 J 
2c 2c 


(17-c) 


A/ _ _ A/ r 

tan0 4 — — (£ 3 + £ 4 ) + tan0 3 — — I £ 0 + 2(£j + £ 0 + £ 3 ) + £ 4 I 
2c 2c 


(17-d) 


A/ _ Air __ _ _ _ 1 

tan 0 n - —(£„_! + £ n )+ tan @ n _\ - — [£ 0 + 2(£, +£ 2 +£ 3 + .... + £ n _\) + £ n \ 
2c 2c 


A / 
2c 


n—l 

£o + 2 X + 

i=l 


(17-e) 


Modified Deflection Equation 

By replacing the bending strain, £ ( . , in deflection equation (12) with the true bending strain, 
£, , one obtains the following modified deflection equation: 


^2 


3^3 = 


(A If 3 

(2£ 0 +£,) 
6c 

(18-a) 

(A If , _ 3 

(5£ 0 + 6f, + £ 2 ) 
6c 

(18-b) 

6c [>£o + 6(2£, + £ 2 ) + £ 3 ] 

(18-c) 
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(18-d) 


^4 = 


(A/) 2 r _ / _ _ _ 

— — — ^1 1£ 0 + 6( 3 £| + 2£ 2 + £3 



(A IT r r _ _ _ , _ 

y n - — |(3n- 1)£ 0 + 6[(«- l)£j +{n- 2)e 2 + (n - 3)£ 2 + .... + [n- (n — l)]£ n _i J + £ n \ 

DC 

„-l _ _] ( 18 ' e ) 

(3 n - 1)£ 0 + 6 ^ (n - /)£,• + £ n 
i=l 

When the true bending strain data {£ 0 ,£ p £ 2 ,£ 3 ,£ 4 ,.. ..,£„} are substituted, the slopes and deflections 
of the tubular spar at the strain-sensing stations can be calculated from equations (17) and (18), 
respectively, for the combined bending and torsion loading case. 

TAPERED CANTILEVER BEAMS 

For a cantilever beam with a depth that tapers down arbitrarily from the beam root to the 
beam tip (for example, a conventional aircraft wing), the half beam depth, c(x) , is no longer 
constant along the span. Assuming that c(x) varies linearly in the region .y. , <x<x { between the 
two adjacent strain-sensing stations, {r H ,r f }, c(x) can be expressed as 


(A/) 2 

6 c 


x — X 

c(x) = c i _ l -(c i _ 1 -c i ) — ; Xj_i < x < Xj (19) 

A l 

in which {c ; _,,c ( } are the values of c(x) at the strain-sensing stations, { , , x ( }, respectively. 

If the beam depth tapers down linearly from the beam root to the beam tip, then { c i _ l ,c i } can be 
expressed as 


, 2 0'-l)A/ , ,(/-!) , \ / 

c i _ l =c 0 -(c 0 -c n ) = c 0 -(c 0 -c n ) ; / = nAl 

/ n 


( 20 ) 


i(Al) i 

Ci=c 0 -(c 0 -c n )—— = c 0 -(c 0 -c n )- ; / = nAl 

l n 


( 21 ) 


in which jc 0 ,c n } are the values of c(x) at the beam root and beam tip, respectively. 

Slope Equations 

The slope, tan 0(x ) , of the tapered beam in the region, x i _ { <x< x t , between the two adjacent 
strain-sensing stations, [x t ,_ p jc ; } , can be obtained by integrating equation (3) (c is replaced with 
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c(x) ), with constant of integration determined by enforcing the continuity of slope at the inboard 
adjacent strain-sensing station, x i _ { , as 


tan0(x) = [ 

J x, 


d 2 y 

— jdx+ tant/,_! 


dx l 


Slope increment 


Slope at Xj_ 


£(x) 

c(x) 


dx + tan0, 


i - 1 


A-i < x < x,. 


( 22 ) 


By substituting strain equation (5) and half-depth equation (19) into slope equation (22), and 
carrying out the integration, one obtains the slope, tan (9 = [tan6>(x, )] , at the strain-sensing station, 
x, , as (see Appendix D for details of integration) (ref. 4) 


tan 6 = A l 


£; 


i - 1 


' _j_ £ i-l C i 


(C.--1-C,-) 


£,C ‘ 1 log — 


+ tan 6 


(i = 1, 2, 3, n) 


(23) 


Deflection Equations 

The deflection, y(x) , of the tapered beam in the region, x ( _, <x<x,, between the two 
adjacent strain-sensing stations, {x ; _,,x,}, can be obtained by integrating slope equation (22), 
with the constant of integration determined by enforcing the continuity of deflection at the inboard 
adjacent strain-sensing station, x ( _, , as 


r x r* c x ^(•^-) r 

y(x)= tan 0(x)c/x + y jl = dxdx + 

J X: 1 l—,, 1 J X; , J X: . r( S' \ * 


Integration of slope 


Deflection 
at x,_i 


c(c) 


tan0, y/x + 


Deflection increment 


Deflection at x 
due to tan0._i 


U-\ 
Deflection 

at X:_t 




(24) 


By substituting strain equation (5) and half-depth equation (19) into deflection equation (24), and 
carrying out the integration, one obtains the deflection, y. = [ y(x ( ) ] , at the strain-sensing station, 
X,. , as (see Appendix D for details of integration) (ref. 4) 


t = ( a /) 2 


£ i - 1 £ i £ i-\ C i £ i C i - 1 


2(c,. 1 -c,) (<:-•,_! - ) 


c, log — + (£)-!-£■, ) 


+ y,_, + A l tariff , 


(25) 


(i = 1,2, 3,....,n);y 0 =tan0 o =0 


SLIGHTLY TAPERED CANTILEVER BEAMS 

For a slightly tapered cantilever beam, the perturbation and stepwise methods can be used to 
develop the slope and deflection equations. This section discusses these two methods. 
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Perturbation Method 


When the cantilever beam is slightly tapered (that is, (c ; /c,._,) — > 1), the logarithmic terms 
in slope equation (23) and deflection equation (25) will approach zero [that is, log(c, /<:,_,) — > 0]. 
Therefore, the perturbation method must be used to expand the logarithmic term, log(c ( . /<?,._,) , 
in the neighborhood of (c i /c i _ 1 ) ~ 1 and then obtain nonzero mathematical expressions so that the 
slope and deflection equations could degenerate into those for the uniform cantilever beam case, 
(c i /c i _ l )=\. 

Slope Equations 

For a slightly tapered beam, the logarithm term, \og(c i /c i _ l )^l, in slope equation (23) 
can be expanded in the neighborhood of (c i /c j _ 1 ) ~ 1 [that is, log(c ; /c I _ 1 ) ~ 0 ]. Carrying out the 
expansion up to the second order terms in(c, t - c, ) 2 , the logarithm term, log (c ; /c ; _, ) , takes on the 
following form as c ; approaches \{c i /c i _ l ) — > l] (see Appendix E for details of mathematical 

expansions) (ref. 4): 



2c 2 


-(c,. - 3c w ) 



(26) 


In light of equation (26), slope equation (23) becomes the following form for the slightly tapered 
beam: 


tan 0 j = 


A/ 

2c, ,i 



+ tan(7 , 



(27) 


When the tapering diminishes (that is, c, = c,_, = c ), equation (27) degenerates into the slope 
equation for the uniform cantilever beams as 

M 

tan^,. =—(e,.-i +£,.) + tan 0 ( ._i ; c, = c,_, = c (28) 


which is exactly the same as slope equation (7) for the uniform beam. 

Deflection Equations 

When the beam becomes slightly tapered (that is, c, — > c i _ ] , the term, 
[(■. logic,. /c, ,) + (c, , in equation (25) can be expanded in the vicinity of (c, / c, ,) = 1 (that 

is, logic /c, , ) = ()). Carrying out the expansion up to the third order terms of (c, , - c, ) 3 , the term, 
[c, log( ,) + (c, , - c,)] , takes on the following form as c ( . approaches c, , [(c,/c, ,) — > l] (see 
Appendix E for details of mathematical expansions) (ref. 4): 


15 



c,-log(— ) + (c i _ 1 -c i ) 

c i - 1 


( C /-1 - c , ) 2 


6 c f _, 


3 + 


i.Ci-1-Ci) 


c , 




->1 


C /-l 


(29) 


In light of equation (29), deflection equation (25) becomes the following form for the slightly 
tapered cantilever beam: 


y> = 


(A If 


6c,_i 


3-— K, +£, 

c /-i 


+ y,_[ + A/tan6f_, 


->1 


(30) 


"/—l 


When the beam taper diminishes (that is, c, = c i _ l = c ), equation (30) degenerates into the deflection 
equation for the uniform cantilever beam as 

(A If 

y, = — — (28^ + £,)+ V / , + A/ tan 0,_, ; c,. = c,_, = c (31) 

6c 

which is identical to the last expression of deflection equation (10) for the uniform cantilever 
beam. 


Stepwise Method 

The slightly tapered cantilever beam can also be approximated with a stepwise tapered 
cantilever beam. The beam segment between the two adjacent strain-sensing stations, 
{x t , x t }, can be considered as a uniform beam segment with an average constant semi-depth 
of c = (1/2 )(c ( ._j + c ; .) . Slope equation (8) and displacement equation (12) can be modified for the 
stepwise tapered cantilever beam as follows. 

Slope Equation 

The slope equation of the stepwise tapered cantilever beam can be written as 


tan 6 { = A / 


gp+g, 

_c 0 +q 


+ tan0 o ; tan0 o =O 


(32-a) 


tan 0 2 = A l 


tan0, = A l 



+ tan = A l 


+ tan 0 2 = A l 


£i g 2 _|_ g« gj 

Ci+C 2 Cp+Cj 


+ tan 6, t 


g 2 + g 3 ! £l + g 2 j g P + g I 

C, + C 3 C[ + C 2 C () + Cj 


I + tan 6„ 


(32-b) 


(32-c) 
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tan 0 4 = A l 


£3 + g 4 

C3 + C4 


+ tan0, = A l 


g 3 + g 4 t g 2 + g 3 t g[ + g 2 | g() + g l 

C3 + c 4 c, + c 3 C| + c 2 c 0 + Cj 


+ tan0 n 


tan0„ = A/ 


g„- 1 + g„ 

c »-i + c » 


+ tan0 ;j _, = A/^ 


g/-i + g, 

c i-i + c ,- 


+ tan0„ 


Deflection Equations 

The deflection equation of the stepwise tapered cantilever beam can be written as 


A/ 2 , , A/ 2 

= — — ; — :(2g 0 + e, ) + .Vo + Al tan ^o = — 
3 (c 0 +q) 3 


2e 0 + £j 

c 0 + C, 


; y 0 =O,tan0 o =0 


A/ 2 , , A/ 2 

* = 37 T 71 (2£ ' +e3+y ' + A/tan0 ' = T 


2 g t + g 2 | 2 e 0 + e. 


c, +c 2 


C 0 +Ci 


+ A l tan0. 


T 3 = 


A/ 2 


3(c 2 + c 3 ) 
+A/(tan0[ + tan0 2 ) 


(2£ 2 + g 3 ) + T2 + Al tan0 2 - — - 


2g 2 + £3 | 2£, + £ 2 

c 2 +c 3 Cj +c 2 


A/ 2 , , A/ 2 

V 4 = — — ; — “(2g, + e 4 ) + y 3 + A/ tan 0 3 = — 


3(c 3 + c 4 ) 

+ A / ( tan 0, + tan 0 2 + tan 0, ) 


2 fi 3 + £ 4 2 £ 2 + £3 

L C 3 +C 4 C 2 +C 3 


A/ 2 , 3 A/ 2 

— : — :(2 g „ 1 + g J + y „- . + A/ tan0 « . = — 
3 ( C «-1 + O 3 


2g«-i + g„ | 2 £„_ 2 + £„_1 

C «-l + C „-2 + C „-l 


+A/^tan 6 I 


The last term in equation (33-e) can be expressed as 


(32-d) 

(32-e) 

(33-a) 

(33-b) 

( 33 -c) 

(33-d) 

( 33 -e) 
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n — 1 


Xtan0,.=A/ ( n -l)^L + ( n _2)^l + (n -3) e2 + e 3 . . r. 1 ^.- 2 +<W 


c 0 +c. 


Cl +c 2 


c 2 +c 3 


- + .... + [« — (n — 1)]- 


C n - 2 + C„_! 


«— 1 


= A/X 


(n - i ) 


c,-_i + c ; 


( 34 ) 


In light of equation (34), deflection equation (33-e) becomes 


y n = ' 


(A If 


2g «- 1 + g « | 2e „- 2 + g «-i 


+ c n 


C n - 2 + C n - 1 


+ 


(A/) 2 X 


(n - i ) 


g /-i + g / 

C/-i + c,- 


(35) 


Note that equation (35) is for the last span-wise strain-sensing station, n. It can, however, 

be used to calculate the deflection at any strain-sensing station, v, (= 1,2,3, , n), by simply 

replacing the indices, {n, /}, with {/, /}. 


OTHER APPLICATIONS 


The bending deflection equations developed for the uniform and tapered cantilever beams 
in the previous sections can also be applied to predict the deformed shapes of other types of 
structures. This section discusses the application of these equations to cantilever wing boxes, 
two-point supported beams (simple beams), and plates. 

Wing Boxes 

Figure 4 shows a wing box, which is a type of a tapered cantilever beam. The simple way 
to monitor the wing box deformation under combined bending and torsion is to use two parallel 
bending strain- sensing lines in the span- wise direction (fig. 4). If { y j , y' } denote the deflections 
of the front and rear sensors, respectively, at the same span-wise strain-sensing cross section, .x / 

(i = 1, 2, 3, n), and if cl is the chord- wise distance between the front and rear sensor lines, then 

the cross-sectional twist angle, 0, , at any strain-sensing cross section, x t (fig. 4), can be calculated 
from the following equation: 


0 ,- 



(36) 


In equation (36), { y t , y' } can be calculated from deflection equation (25) (or eq. (30) or (35)) 
for a tapered cantilever beam by using the strain data obtained from the two-line sensing system as 
inputs. With equation (36) on hand, the distortion sensor is not required. 
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Two-Point Supported Beams 


Deflection equation (12) was developed for a uniform cantilever beam with a clamped left 
end (that is, zero slope and zero deflection, y n = tan 9 0 = 0 ) and a finite deflection ( y n ) at the beam 
tip. If the strains, y n = 0 , of the two-point supported beam with a clamped left end are used as 
inputs, then equation (12) theoretically should give zero deflection (y n = 0 ) at the right end. If the 
calculated right end deflection is not zero, however, and it has a small value because of possible 
measurement errors, then equation (12) must be corrected. 


To eliminate the nonzero (y n ^ 0 ) right end deflection, the deflection curve is rotated slightly 
as a whole, with respect to the left end support point, in such a way as to bring the tip deflection, 
y n (eq. (12-e)), down to zero ( y n = 0 ). Therefore, the deflection calculated from equation (12) for 

each strain-sensing station, x i , is reduced by an amount of ( i/n)y n (i = 1, 2, 3 n), a linearly 

varying fraction of the nonzero right end cantilever tip deflection, y n . 


If yf denotes the deflection of the two-point supported beam at the strain-sensing station, v, , 
then the deflection equation for the two-point supported beam can be established from cantilever 
deflection equation (12) with corrections as follows: 


1 


Ti 


yi~-y n 
(A/) 


n 

\2 


1 


(2e 0 + e l ) — y n 
6 c n 


(37-a) 


^ =yi~-y n 

n 


(A/) 2 , v 2 

— — (5e 0 + 6e, +e 2 ) — y n 
6c n 


y 3 =y 3 --y n 

n 


(A/) 


[8e 0 +6(2e, +e 2 ) + e 3 ] — y n 

6C 71 


3^4 = 3 7 4 y „ 

n 


~ ~ [ILo + 6(36, +2e 2 +e,) + e 4 ] — y n 
be n 


(37-b) 


(37-c) 


(37-d) 
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(37-e) 


B n 

y„ =y n --y n 

n 

(A/) 2 T vC \ 1 n 

= —— ('in-l)e 0 +6Y J {n-i)e i +e n -~y„ = 0 
6 c l ;=i J n 

In equation (37), y t (i = 1, 2, 3 n) represents the deflections calculated from cantilever beam 

deflection equation (12). 

If the two-point supported beam is symmetrically loaded, the distribution of the strain, e, , 
will be symmetrical with respect to the point of symmetry, and the mathematical expressions of 
the deflections obtained from equation (37) can be proven to be identical for the two mirror image 
sensing stations. Deflection equation (37) established for the two-point supported beam can also 
be used to predict the deformed shapes of plates. 


Beam With Two Ends Clamped 

For the beam with two ends clamped, the support conditions for the two ends are as follows: 


At the left clamped end: 

=0 

o 

II 

<3? 

§ 

(38) 

At the right clamped end: 

y„ = o ; 

tan 0 n = 0 

(39) 


When the strain values, £, , for the present case are used, beam deflection equation (37) should 
produce condition (39) (y„ = tan 0 n = 0 ) at the right clamped end. 

Clamped and Simply Supported Beam 

The two-point supported beam with the left end clamped and right end simply supported is a 
special case of a cantilever beam with the tip deflection, y n (fig. 2) (eq. (12-e)), forced to zero. The 
end support conditions for the clamped and simply supported beam are as follows: 


At the clamped end: 

=° 

; tan0 o = 0 

(40) 

At the simply supported end: 

J „=0 : 

; tan0 (l ^ 0 

(41) 


When the strain values, £ ( , for the present case are used as inputs, beam deflection equation (37) 
results in condition ( 41 ) ( y n ~ 0 , tan 0 n y- 0 ) at the right simply supported end. 

Beam With Two Ends Simply Supported 

For the beam with two ends simply supported, the beam deflection slope at the simply 
supported left end is no longer zero ( tan 6 0 ^ 0 ). Therefore, to formulate the slope and deflection 
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equations for this beam, the nonzero term, tan 0 0 0 (which was neglected for the cantilever beam 

case), must be retained in both slope equation (8) and deflection equation (12). 

Slope Equation 

The slope equation for the beam with two ends simply supported can be written as 
Al 

tan 6>, = — (e 0 + £, ) + tan 0 0 ; ( tan 6 0 ^ 0 at left support) (42-a) 


tan 0 2 


(£| + £ 2 ) T tan 

2c 


A / 
2c 


(£ 0 + 2£, + £ 2 ) + tan0 o 


(42-b) 


tan£h 


A l_ 
2c 


A / 


(£., + £j) + tan Or, — — ^£ 0 + 2(£, + £ 2 ) + £3 j + tan 0 0 


2c 


(42-c) 


A / Air n 

tan0 4 — — (£3 + £ 4 ) + tan0 3 — — I £ 0 + 2 (£j + £ 2 + £3) + £ 4 I + tan0 o 
2c 2c 


(42-d) 


A / Air 1 

tan0 « - ^(£,,-1 + e «) + tan0 n _, - — [e 0 + 2 ( e i + e 2 + e 3 + .... + £„_!) + £„] + tan 0 O 


A/ 

2c 


n — 1 


(42-e) 


e o + 2 X e / + e » 


+ tan0„ 


(=1 


Deflection Equations 

Let yf be the deflection of the two-point supported beam at the strain-sensing station, x t . 
Deflection equation (12) can then be modified for the beam with two ends simply supported as 


yf = y, + A l tan 0 Q 

AM 1 ( 2 £ „ + £| ) + Aaa„e 0 ; ( - v o = e o = 0. tan^ ^ 0 at left support) 


(43-a) 


yf = y 2 + 2 Al tan 0 0 
(A/) 2 


6c 


(5£ 0 + 6 £j + £ 0 ) + 2 Al tan 0 O 


(43-b) 


yf = y 3 + 3 Al tan 0 0 


(A/) 2 

6c 


[8£ 0 + 6(2£, + £ 2 ) + £ 3 J + 3 Al tan0 o 


(43 -c) 
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(43-d) 


yl = y 4 + 4 A / tan0 o 


(A/) 2 r / 

— — — [^1 l£ 0 + 6 (3c, + 2e 2 + £ 3 


) + £ 4 J + 4 A l tan 0 0 


y B n = y„ + nA 1 tan 6 0 


(A/) 2 

6c 


n— 1 


(3n - 1)£ 0 + 6^(77 - i)£. + £„ 


+ 77 A/ tan O n 


(£0 = £ n = 0) 


(43-e) 


in which y ( . (7 = 1, 2, 3 n) is the deflection calculated from cantilever beam deflection 

equation ( 1 2). By setting y K n = 0 in equation (43-e) for the right simply supported end, the unknown 
slope, tan 0 0 , at the left simply supported end can be expressed in terms of the measured strains, £, , as 


tan0 o 


1 (A If 
n(Al) 6c 


n - 1 

(3/7 - 1)£ 0 + 6^(77 - i)e. + e n 

7 = 1 


1 

n(Al) yn 


(44) 


In equation (44), the last term (containing y n ) can be obtained from equation (12-e). In light 
of equation (44), beam deflection equation (43) takes on the form of equation (37). Thus, for 
the beam with two ends simply supported, equations (43) and (37) will result in identical beam 
deflection curves. 


Symmetrically Loaded Beam 

If the beam is symmetrically loaded (concentrated or distributed), then tan 0 O = - tan 0 n , 
and equation (42-e) becomes 


A/ r 1 

tan0 o -~~[ e o +2 (£i +£ 2 +e 3 + .... + £„_,) + £„] 


( £ 0 


0 ) 


(45) 


For a beam that is symmetrically loaded and simply supported, the slopes, tan 0 () , calculated 
from equations (44) and (45) should be identical, although their mathematical expressions are 
different. This match can be verified by considering an example of a simply supported beam 
in which n = 4 ( e t] = £ 4 = 0 ) under central vertical load. The symmetry condition results in 
( £3 = £j ), and equation (44) becomes 


tan 0 O 


— [1 1 £q + 6(3£, + 2fi 2 + £3) + £4] 

24c 

A l 

— ( 2 £, +£ 2 ) 

2 c 


(n = 4; £ 0 = £ 4 = 0 ; £3 = £j) 


(46) 


When the same strain conditions are applied ( £ 0 = £ 4 = 0; £ 3 = £, ), slope equation (45) becomes 
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; (n — 4; e 0 — e 4 — 0 ; e 3 — e, ) 


(47) 


Mr n 

tan 0 0 — — — — ^£ 0 + 2(£| + £, + £j ) + £ 4 J 
M 

- “~( 2e i + e 2) 

2c 

For the present centrally loaded simply supported beam (n = 4; £ 0 = £ 4 = 0 ; £ 3 = £, ), 
equations (46) and (47) result in identical mathematical expressions for the left end slope, 
tan 9 0 ■ 


Plates 

For plates, the strain-sensing system requires multiple parallel strain-sensing lines (fig. 5) 
across the two opposite edges to generate a three-dimensional deformed shape. A narrow strip 
of the plate along each sensing line can be considered as a two-point supported beam. Hence, 
deflection equations (37) or (43) developed for the two-point supported beam also are applicable 
to the predictions of the deformed shapes of plates. 

VALIDATION OF THEORY 

This section verifies the mathematical accuracy of the bending and torsion displacement 
equations developed for the uniform cantilever beam. The mathematical expressions of the 
displacement equations are compared with those of classical cantilever beam theory. 

Bending Case 

To test the accuracy of the mathematical expressions of the slope and deflection equations, 
a cantilever beam of length l subjected to an upward tip load of P is used as an example. For this 
case the moment, M, and associated bending strain are linearly decreasing functions of x from the 
built-in end of the beam toward zero at the beam tip. 

Assume that the cantilever beam has five bending strain sensors { £ 0 ,£j,£ 2 ,£ 3 ,£ 4 } equally 
spaced over the span by M = 1/4 along the bottom surface of the beam. The purpose is to determine 
whether the calculated functional expressions of the tip slope and tip deflection of the cantilever 
beam match those given by classical beam theory. 

Bending Strains 

For a cantilever beam of length l subjected to a tip load of P, the bending strain (or bending 
moment) is a linearly decreasing function of x from the built-in end of the beam down to zero at 
the free end of the beam. The bending strains at the strain-sensing stations, x t (=0, 1, 2, 3, 4), can 
then be written as 
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Tip Slope 

When strain equation (48) is substituted into slope equation (8-d), the beam tip slope is 
calculated as 


tan 0 = tan 0, 
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M(0)c 

in which the moment-strain relationship, e o = (eq. (2)), withM(O) = PI at the built-in end 

El 

( x = 0) is used. Equation (49) is identical to the classical tip slope expressions for a cantilever beam 
of length l subjected to a vertical load of P at the free end (refs. 2, 3). 

Tip Deflection 

When strain equation (48) is substituted into deflection equation ( 1 2-d), the beam tip deflection 
is calculated as 


24 



y « = 3^4 = 


(A/) 2 r / 

1 1£ 0 + 6 ( 3 £| +2 £ 2 + £ 3 

6c 



fo. 

6c 


11 + 6 


_ 3 _ 2 n 

3x— +2x— +— +0 
4 4 4 J 


_ / 2 M(0)c _ l 2 ( Pl)c 
3c 0 _ 3c El ~ 3c El 
PI 3 
3 El 


( 50 ) 


M(0)c 

Again, the relationship g = — — — (eq. (2)) at the built-in end (x = 0) and M( 0) = PI is used. 

El 

Equation (50) is identical to the tip deflection expressions for a classical cantilever beam of length 
l subjected to a vertical load of P at the free end (refs. 2, 3). 

The present example verifies the accuracy and dependability of the slope and deflection 
equations that have been developed. This example also provides confidence for applying these 
displacement equations to monitor the in-flight deformed shape of the Helios flying wing. 

Torsion Case 

Assume that the cantilever beam has four distortion strain sensors { 7 0 , 7, , / 2 , 7.3 } 
equally spaced by Al = 1/4 along the cantilever beam, which is subjected to a twisting moment 
of T at the tip. The purpose is to determine whether the calculated functional expression of the tip 
twist angle matches that given by classical beam theory. 

Distortion Strains 

For the case of constant twisting along the span (that is, T 0 = 7j = T 2 = Tj = T 4 = T ), the 
surface twist angles (shear strains) are constant at all sensing stations. Namely, 


7i = 7 2 = 7 3 = 7 o 


(51) 


Tip Twist Angle 

From equation (15-d), the tip twist angle, 0 n (= 0 4 ) , can be written as 

A / 

t = 04 = — (7o + 7 1 + 7 2 + 7 3 ) (52) 

c 

In light of equation (51), the tip twist angle of equation (52) becomes 
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( 53 ) 
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0„ = 04 = — (r 0 + 7i + r 2 + r 3 ) = — (4y 0 ) = -y 0 
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Equation (53) is exactly the same as the expression for the tip cross-sectional twist angle of the 
classical cantilever beam of length l subjected to a tip twisting moment of T (refs. 2, 3). 

FINITE-ELEMENT ANALYSIS 

To test the prediction accuracy of deflection equation (12), strain sensor data and measured 
deflection data are needed. Before the experimental data are available, an alternative approach 
is to use finite-element analysis to generate the desired bending strain and deflection data. The 
Structural Performance And Resizing (SPAR) finite-element computer program (ref. 5) has been 
used for this purpose. The next section describes the method for generating the bending strains, e, , 
from the SPAR outputs. 


l_( 7b 
c{gJ 


Generation of the Bending Strains 

From the SPAR nodal displacements outputs, the axial displacements, {0 P 0 2 } , at the two 
ends of the element lying in the lower outermost surface can be used to calculate the bending strain 
for that particular element. If the strain-sensing station, x t , is located in the middle of the element, 
then the displacement differential, (<5 2 - <5, ) , can be divided by the finite-element span-wise length, 
e, to generate the desired bending strain, £, , at the strain-sensing station, x t . Namely, 


&-g, j 


(54) 


If the strain-sensing station, x t , is located at the juncture of two adjacent elements, then the 
following average strain value is used for the strain, £ i , at the strain-sensing station, x t : 
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8 1 -8 i 
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&-g, j 


(55) 


In equation (55), { ( ) i+ , ( ),_ } denote the elements on the outboard and inboard side, respectively, 
of the strain-sensing station, x t . 

The SPAR bending strains calculated from equations (54) and (55) can be input to 
equation (12) to calculate the beam deflections. These calculations can then be compared to the 
deflections obtained from the SPAR outputs to check the prediction accuracy of displacement 
equation (12). 
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Cantilever Tube Models 


Figure 6 shows a typical finite-element model for a cantilever tube with a length of 
l =100.5 in. The model consists of 100 four- node shell elements in the axial direction and 36 in the 
circumferential direction. The end disk was added so that the tip load can be applied at the disk 
center. The cantilever tube is fixed in space at the left end and subjected to an upward vertical load 
of P = 100 lb at the tip. The following three types of tubes fabricated with different materials were 
considered: 

1. Aluminum tube 

2. Four-ply composite tube: (45/-45/-45/45) carbon fiber composite tube 

3. Helios composite tube: 4 plies (45/-45/-45/45) of carbon fiber composite wall with 
27 plies (0 8 /90/0 g /90/0 9 ) of carbon fiber composite reinforcements in the 60-deg regions 
of the top and bottom of the spar caps (actual Helios wing spar). 

Table 1 lists the dimensions and table 2 lists the material properties of these types of tubes. 


Table 1. Dimensions of three types of tubular spars. 


Tube type 

/, in. 

a, in. 

c, in. 

t, in. 

Aluminum 

100.5 

4.0 

4.01148 

0.02296 

Four-ply composite 

100.5 

4.0 

4.01148 

0.02296 

Helios composite 

100.5 

4.0* (4.07749)** 

4.16646** 

0.02296* (0.17794)** 


* Four-ply regions, ** 31 -ply regions 


Table 2. Material properties of carbon fiber composite and aluminum tubes. 


Carbon fiber composite 

Aluminum 

E l = 16.3 x 10 6 lb/in 2 

E = 10.5 x 10 6 lb/in 2 

E t = 0.98 x 10 6 lb/in 2 

NA 

G lt = 0.71 lb/in 2 

G = 4.0 x 10 6 lb/in 2 

v LT = 0.34 

v = 0.33 


Wing Box Model 

Figure 7 shows the tapered cantilever wing box model. The box is made of aluminum 
(table 2), and has a length of l = 100 in., width of w = 17 in., and depth of 2c 0 =4 in. at the 
root, which tapers down to a depth of 2 c n =3 in. at the tip. The upper and lower skins and 
the vertical walls have an identical thickness of t = 0.02296 in. (wall thickness of an aluminum 
cantilever tube). The wing box was subjected to a vertical load of P = 100 lb and a torque of 
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T = (Pi 2) x w = 850 in-lb at the wing box tip. The torque was generated by applying two equal and 
opposite forces of Pi 2 = 50 lb at the two comers of the wing box tip (fig. 7). 

Beam Models 

The two-point supported beam models are identical to the cantilever tube model (fig. 6) 
except that the support conditions have been changed. Figure 8 shows the two-point supported 
tubular beam models subjected to a downward force of P = -200 lb at the beam center, under three 
cases of support conditions: (1) both ends clamped (fig. 8(a)), (2) left end clamped and right end 
simply supported (fig. 8(b)), and (3) both ends simply supported (fig. 8(c)). The deformed shapes 
also are shown. 


Plate Model 

Figure 9 shows the dimensions of an aluminum square plate (/ = 24 in., w = 24 in., 
t = 0.2 in.). The plate was subjected to a downward force of P = -400 lb at the plate center, 
under two edge-support conditions: (1) four edges clamped (fig. 9(a)), and (2) four edges simply 
supported (fig. 9 (b)). Because of loading symmetry, only one quarter of the plate was modeled. 
For the quarter model, a downward force of only P/4 = -100 lb was applied. 

RESULTS 

The SPAR displacement outputs were used to generate the bending strain, e l , for each strain- 
sensing station. These strain values were then input to the appropriate deflection equations to 
calculate the deflections at all the strain-sensing stations of all the study cases. 

Cantilever Tubes 

This section presents the results of three types of cantilever tubes (aluminum, four-ply 
composite, and Helios composite). The deflections calculated from the Ko displacement theories 
are compared with those calculated from the SPAR finite-element analysis. 

Aluminum Cantilever Tube 

Figure 10 shows the bending strains, e, , at the strain-sensing stations on the aluminum tube, 
calculated from the SPAR nodal displacement outputs. The strain data from figure 10 were used to 
calculate the deflections, y t , from deflection equation (12). To find the minimum number of strain 
sensors required for acceptable accuracy of the deflection predictions, three cases of n (2, 4, 8), 
were considered. Table 3 compares the deflections, y t , calculated from deflection equation (12) to 
those obtained from the SPAR outputs for the different n values. Note that the difference between 
the deflection values calculated from the SPAR program and those calculated from deflection 
equation (12) is minimal. Also note that the effect of the number of sensing stations, n, is negligible 
for the present uniform cantilever aluminum tube. 
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Table 3(a). Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (12); aluminum cantilever tube; n = 2. 


Deflection ( y t ), in. 


% 

3b 

3b 

SPAR 

0.0000 

0.2253 

0.7138 

Deflection equation (12) 

0.0000 

0.2184 

0.7000 

Difference, percent 

0.0000 

3.0451 

1.9389 


Table 3(b). Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (12); aluminum cantilever tube; n = 4. 


Deflection ( y t ), in. 


% 

A 

33 

3b 

34 

SPAR 

0.0000 

0.0636 

0.2253 

0.4530 

0.7138 

Deflection equation (12) 

0.0000 

0.0604 

0.2206 

0.4468 

0.7056 

Difference, percent 

0.0000 

5.1242 

2.0996 

1.3707 

1.1544 


Table 3(c). Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (12); aluminum cantilever tube; n = 8. 





Deflection ( y t ), 

in. 






3b 

3 7 i 

3b 

3b 

3b 

3b 

3b 

3b 

3b 

SPAR 

0.0000 

0.0179 

0.0636 

0.1341 

0.2253 

0.3330 

0.4530 

0.5814 

0.7138 

Deflection 
equation (12) 

0.0000 

0.0157 

0.0605 

0.1302 

0.2205 

0.3273 

0.4465 

0.5738 

0.7051 

Difference, 

percent 

0.0000 

12.2142 

4.9356 

2.9223 

2.1262 

1.7208 

1.4524 

1.3072 

1.2230 


Data from table 3 are plotted in figure 11 for different n values. The two sets of deflection 
curves are quite close, indicating the accuracy of displacement equation (12). The deflection curves 
calculated from deflection equation (12) lie slightly below those obtained from the SPAR outputs. 
Notice that the difference in the tip deflection did not improve when n was increased from 4 to 8. 
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Four- Ply Composite Cantilever Tube 


Figure 1 2 shows the bending strains, e j , at the strain-sensing stations on the four-ply composite 
tube calculated from the SPAR outputs. In this case the SPAR strain curve is practically linear. The 
strain data from figure 12 were used to calculate the deflections, y t , from deflection equation (12) 
under three cases of n (2, 4, 8). Table 4 compares the deflections, y n , calculated from deflection 
equation (12) to those obtained from the SPAR outputs. Again, the difference between the 
values calculated from the SPAR program and those calculated from deflection equation (12) 
is minimal. 


Table 4(a). Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (12); four-ply composite cantilever tube; n = 2. 


Deflection ( y t ), in. 


y 0 

A 


SPAR 

0.0000 

0.9110 

2.8946 

Deflection equation (12) 

0.0000 

0.9170 

2.9373 

Difference, percent 

0.0000 

0.6641 

1.4741 


Table 4(b). Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (12); four-ply composite cantilever tube; n = 4. 


Deflection ( y t ), in. 


% 

A 

>2 

y 3 

y 4 

SPAR 

0.0000 

0.2577 

0.9110 

1.8347 

2.8946 

Deflection equation (12) 

0.0000 

0.2515 

0.9156 

1.8553 

2.9331 

Difference, percent 

0.0000 

2.4058 

0.5116 

1.1223 

1.3290 
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Table 4(c). Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (12); four-ply composite cantilever tube; n = 8. 


Deflection ( y t ), in. 


% 

34 

35 

34 34 

34 

34 

34 

34 

SPAR 

0.0000 

0.0749 

0.2577 

0.5421 0.9110 

1.3475 

1.8347 

2.3558 

2.8946 

Deflection 
equation (12) 

0.0000 

0.0658 

0.2519 

0.5413 0.9167 

1.3610 

1.8570 

2.3877 

2.9356 

Difference, 

percent 

0.0000 

12.2264 

2.2545 

0.1531 0.6301 

1.0000 

1.2176 

1.3526 

1.4154 


Data from table 4 are plotted in figure 13 for different n values. The two sets of deflection 
curves are quite close, indicating the accuracy of deflection equation (12). The deflection curves 
calculated from deflection equation (12) lie slightly above those obtained from the SPAR outputs. 
Again, the effect of n is negligible. 

Helios Composite Cantilever Tube 

Figure 14 shows the bending strains, £, , at the strain-sensing stations on the Helios composite 
tube calculated from the SPAR outputs. In this case the SPAR strain curve is nearly linear. The 
strain data from figure 14 were used to calculate the deflections, y t , from deflection equation (12) 
under three cases of n (2, 4, 8). Table 5 compares the deflections, y l , calculated from deflection 
equation (12) to those obtained from the SPAR outputs. Again, the difference between the 
values calculated from the SPAR program and those calculated from deflection equation (12) is 
minimal. 


Table 5(a). Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (12); Helios composite cantilever tube; n = 2. 


Deflection ( y t ), in. 


34 

34 

34 

SPAR 

0.0000 

0.2167 

0.6721 

Deflection equation (12) 

0.0000 

0.2115 

0.6739 

Difference, percent 

0.0000 

2.3863 

0.2678 
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Table 5(b). Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (12); Helios composite cantilever tube; n = 4. 


Deflection ( y t ), in. 


% 

-V. 

3b 

3b 

34 

SPAR 

0.0000 

0.0673 

0.2167 

0.4284 

0.6721 

Deflection equation (12) 

0.0000 

0.0580 

0.2096 

0.4232 

0.6677 

Difference, percent 

0.0000 

13.7952 

3.2726 

1.1976 

0.6443 


Table 5(c). Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (12); Helios composite cantilever tube; n = 8. 





Deflection ( y ( . ), 

in. 






3b 

3 7 i 

3b 

3b 

3b 

3b 

3b 

3b 

3b 

SPAR 

0.0000 

0.0240 

0.0673 

0.1324 

0.2167 

0.3166 

0.4284 

0.5479 

0.6721 

Deflection 
equation (12) 

0.0000 

0.0152 

0.0577 

0.1235 

0.2089 

0.3099 

0.4227 

0.5432 

0.6676 

Difference, 

percent 

0.0000 

36.7126 

14.2263 

6.6853 

3.5818 

2.1381 

1.3307 

0.8669 

0.6725 


Data from table 5 are plotted in figure 15 for different n values. The two sets of deflection 
curves are quite close, indicating the accuracy of deflection equation (12). The deflection curves 
calculated from deflection equation (12) lie slightly below those obtained from the SPAR outputs 
(especially at lower values of n). Again, the effect of n is insignificant. 

Wing Box 

Figure 16 shows plots of the front and rear bending strains, {£,.,£'}, for the tapered wing 
box obtained from the SPAR outputs. Because of the taper and torsion effect, the two SPAR strain 
curves are no longer linear. The front strain curve is convex downward from the wing root toward 
zero at the wingtip. Conversely, the rear strain curve is convex downward near the wing root and 
then turns to concave downward toward zero at the wingtip. The strain data from figure 16 were 
used to calculate the front and rear deflections, { y t y ' }, from deflection equations (30) (tapered) 
and (35) (stepwise) for n = 8. Tables 6 and 7 compare the front and rear deflections, { y. y ' }, 
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respectively, calculated from equations (30) and (35) ({ n , /} replaced with {/,/}) to those obtained 
from the SPAR outputs. 


Table 6. Comparison of front deflections calculated from SPAR with those calculated from 
deflection equations (30) (tapered) and (35) (stepwise); aluminum wing box; n = 8; front sensors. 





Deflection ( y t ), 

in. 






34 

A 

33 

34 

34 

34 

34 

34 

34 

SPAR 

0.0000 

0.0229 

0.0864 

0.1892 

0.3285 

0.5008 

0.7008 

0.9218 

1.1559 

Tapered 

deflection 
equation (30) 

0.0000 

0.0208 

0.0818 

0.1807 

0.3138 

0.4771 

0.6656 

0.8725 

1.0894 

Difference 

from SPAR, 
percent 

0.0000 

9.3961 

5.3910 

4.5190 

4.4956 

4.7332 

5.0255 

5.3481 

5.7489 

Stepwise 

deflection 
equation (35) 

0.0000 

0.0209 

0.0820 

0.1809 

0.3141 

0.4774 

0.6658 

0.8725 

1.0892 

Difference, 

from SPAR, 
percent 

0.0000 

8.9601 

5.1596 

4.3816 

4.4074 

4.6813 

4.9984 

5.3470 

5.7714 
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Table 7. Comparison of rear deflections calculated from SPAR with those calculated from deflection 
equations (30) (tapered) and (35) (stepwise); aluminum wing box; n = 8; rear sensors. 





Deflection ( y ' ), in. 






/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 


34 

A 

35 

34 

34 

34 

34 

34 

34 

SPAR 

0.0000 

0.0227 

0.0844 

0.1801 

0.3045 

0.4525 

0.6182 

0.7960 

0.9813 

Tapered 

deflection 
equation (30) 

0.0000 

0.0200 

0.0779 

0.1688 

0.2868 

0.4259 

0.5805 

0.7450 

0.9140 

Difference 

from 

SPAR, 

percent 

0.0000 

11.9683 

7.6339 

6.2398 

5.8383 

5.8830 

6.0998 

6.4161 

6.8661 

Stepwise 

deflection 
equation (35) 

0.0000 

0.0201 

0.0781 

0.1690 

0.2869 

0.4259 

0.5804 

0.7446 

0.9133 

Difference 

from 

SPAR, 

percent 

0.0000 

11.5274 

7.4206 

6.1342 

5.7891 

5.8764 

6.1224 

6.4651 

6.9385 


Note from tables 6 and 7 that the difference between the values of { y t y' } calculated from the 
SPAR program and those calculated from equations (30) and (35) is minimal. Tapered deflection 
equation (30) and stepwise deflection equation (35) resulted in extremely close deflections, which 
implies that the tapered and stepwise theories are quite accurate for the present slightly tapered 
cantilever wing box. 

Figure 17 compares the deflections, y, (table 6) and y' (table 7), calculated by means of 
different methods. Note that the deflection curves generated from equations (30) and (35) for each 
sensing line are pictorially indistinguishable and reasonably close to the corresponding SPAR 
deflection curves. 

Table 8 shows the cross-sectional twist angle, 0, , calculated from equation (36) using the 
data of { y t y' } provided in tables 6 and 7, respectively. The data in table 8 show that both tapered 
deflection equation (30) and stepwise deflection equation (35) resulted in extremely close values 
of , which are very close to the values of (f> i calculated from the SPAR program. The similarity 
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among these values implies that either of the two tapered beam theories can be used to obtain quite 
accurate results for the present slightly tapered cantilever wing box. 


Table 8. Comparison of cross-sectional twist angles calculated from SPAR and equation (36) 
using tapered deflection equation (30) and stepwise deflection equation (35); aluminum wing box; 
n = 8. 


Cross-sectional twist angle ( 0,- ), deg 


0o 

01 

02 

03 

04 

05 

06 

07 

00 

SPAR 

0.0000 

0.0008 

0.0070 

0.0308 

0.0809 

0.1629 

0.2784 

0.4239 

0.5883 

Tapered 
deflection 
equation (30) 

0.0000 

0.0027 

0.0130 

0.0399 

0.0910 

0.1727 

0.2868 

0.4299 

0.5914 

Difference 
from SPAR, 
percent 

0.0000 

224.0000 

85.5769 

29.3603 

12.5417 

6.0308 

3.0142 

1.4113 

0.5328 

Stepwise 
deflection 
equation (35) 

0.0000 

0.0027 

0.0131 

0.0401 

0.0915 

0.1735 

0.2879 

0.4312 

0.5929 

Difference 
from SPAR, 
percent 

0.0000 

224.0000 

86.5385 

30.1258 

13.1250 

6.5067 

3.4136 

1.7293 

0.7906 


Note that the high percentage of errors shown in {</>, , </> 2 , </>, } is not significant, because the 
small numbers are divided by small numbers (that is, noise level). The low error levels at the 
wingtip, however, are significant. 

Figure 18 compares the cross-sectional twist angles, 0,- , using the data from table 8. The 
cross-sectional twist angle curves generated by equation (36) using the deflections calculated from 
tapered deflection equation (30) and stepwise deflection equation (35) are pictorially on top of each 
other and very close to the cross-sectional twist angle curve generated from the SPAR program. 
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These results demonstrate that the two-line strain-sensing system could produce quite accurate 
bending deflection and twist angle curves for the present slightly tapered wing box without the 
need for shear strain sensors. 


Simple Beams 

This section presents the results of three types of two-point supported beams (both ends 
clamped, left end clamped and right end simply supported, and both ends simply supported). The 
deflections calculated from the Ko displacement theories are compared with those calculated from 
the SPAR finite-element analysis. 

Both Ends Clamped 

Figure 19 shows the bending strains, £, , for the beam with both ends clamped, calculated 
from equations (54) and (55) using the SPAR displacement outputs. The beam deflections, 
y t , were then calculated from deflection equation (12) using the SPAR strain values presented in 
figure 19. Table 9 presents these deflections and compares them to the deflections, y t , obtained 
from the SPAR outputs and those calculated from classical beam theory. 
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Table 9. Comparison of deflections calculated from deflection equation (12) with those calculated 
from SPAR and classical beam theory; two-point supported aluminum tubular beam with both 
ends clamped; P = -200 lb at beam center; n = 8. 


Deflection ( v. ), in. 


To 

Ti 

t 2 

t 3 

t 4 

(load point) 

T 5 

t 6 

Tv 

Ts 

SPAR 

0.0000 

-0.0047 

-0.0132 

-0.0217 

-0.0264 

-0.0217 

-0.0132 

-0.0047 

0.0000 

Classical beam 
theory 

0.0000 

-0.0034 

-0.0109 

-0.0184 

-0.0218 

-0.0184 

-0.0109 

-0.0034 

0.0000 

Deflection 
equation (12) 

0.0000 

-0.0032 

-0.0107 

-0.0182 

-0.0215 

-0.0182 

-0.0107 

-0.0032 

0.0000 

Difference 
from SPAR, 
percent 

0.0000 

31.9149 

18.9394 

16.1290 

18.5606 

16.1290 

18.9394 

31.9149 

0.0000 

Difference 
from classical 
beam theory, 
percent 

0.0000 

5.8824 

1.8349 

1.0870 

1.3761 

1.0870 

1.8349 

5.8824 

0.0000 


(-) implies downward deflections. 


The deflections, y i , presented in table 9 are plotted in figure 20 for visual comparison. Note 
that the deflection curves generated from equation (12) and those generated from classical beam 
theory agree very well with only a 1.08-percent difference at the beam center. The difference 
between the deflections calculated from equation (12) and those calculated from the SPAR program 
(16.51 percent at the beam center) is attributed to the shear effect, which tends to increase the 
beam deflections. To prove this point, the SPAR program was run again under two cases of input 
conditions. In case 1, the nodal axial displacements were set to zero to eliminate beam cross- 
sectional rotations and obtain deflections purely caused by shear deformation. In case 2, the shear 
modulus, G, was set to infinity (G — >°o) in the SPAR program to simulate the classical beam case. 
Table 10 compares the beam center deflections calculated from the SPAR program under different 
input conditions to those calculated from classical beam theory. 
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Table 10. Comparison of beam center deflections calculated from SPAR with those calculated 
from classical beam theory; two-point supported aluminum tubular beam with both ends clamped; 
P = -200 lb at beam center. 



SPAR 

SPAR 

SPAR 

Classical 


(G*°o) 

(shear only) 

(G— >°o) 

beam theory 

Beam center deflection, in. 

0.2636-1 

0.0442-1 

0.2194-1 

0.2181-1 


The data presented in table 10 prove that the shear effect is indeed considerable for the present 
tubular beam with both ends clamped. In general, the shear effect becomes important only when 
the beam is short with an aspect ratio of lHa less than 4 (lHa < 4, ref. 3, p. 201). The present 
tubular beam has an aspect ratio of lHa = 12.56 (=100.5 -s- 8) and therefore does not belong in the 
short beam category. The thin walled tubular geometry, however, and the two inflection points 
(curvature sign changes) in the deflection curve that divide the beam into three deformed segments 
(fig. 20) seem to cause the beam to exhibit considerable shear effect of approximately 17 percent 
of the total deflection at the beam center. 

Clamped and Simply Supported 

Figure 21 shows the strains, e . , for the beam with a clamped left end and simply supported 
right end calculated from equations (54) and (55) using the SPAR displacement outputs. The 
SPAR strain values presented in figure 21 were used to calculate the beam deflections, y ( , from 
deflection equation (12). Table 11 presents these deflections and compares them to the deflections, 
y t , calculated from the SPAR program and those calculated from classical beam theory. 
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Table 11. Comparison of deflections calculated from deflection equation (12) with those calculated 
from SPAR and classical beam theory; two-point supported aluminum tubular beam with a clamped 
left end and simply supported right end; P = -200 lb at beam center; n = 8. 

Deflection ( » ), in. 



% 

3) 

34 

34 

34 

(load point) 

y 5 

34 

34 

34 

SPAR 

0.0000 

-0.0069 

-0.0200 

-0.0340 

-0.0439 

-0.0415 

-0.0320 

-0.0173 

0.0000 

Classical beam 
theory 

0.0000 

-0.0052 

-0.0170 

-0.0299 

-0.0382 

-0.0376 

-0.0293 

-0.0159 

0.0000 

Deflection 
equation (12) 

0.0000 

-0.0053 

-0.0172 

-0.0300 

-0.0383 

-0.0378 

-0.0297 

-0.0162 

0.0000 

Difference 
from SPAR, 
percent 

Difference 

0.0000 

23.1884 

14.000 

11.7647 

12.7563 

8.9157 

7.1875 

6.3584 

0.0000 

from classical 

0.0000 

1.9231 

1.1765 

0.3344 

0.2618 

0.5319 

1.3652 

1.8868 

0.0000 


beam theory, 
percent 


(-) implies downward deflections. 


The deflection data listed in table 1 1 are plotted in figure 22 for graphical comparison. As in 
the previous case, the deflection curve generated from equation (12) agrees fairly well with that 
generated from classical beam theory with only a 0.36-percent difference at the beam center. It 
is off by 12.79 percent, however, from the SPAR deflection curve at the beam center. Again, the 
difference between the deflections calculated from the SPAR program and those calculated from 
deflection equation (12) is caused by shear deformation. When the shear modulus, G, was set to 
infinity (G — >°°) in the SPAR program, the deflection curve closely approached that calculated 
from classical beam theory. Again, the present tubular beam exhibited the short beam behavior 
(l/2a < 4, ref. 3, p. 201) caused by the existence of one inflection point in the deflection curve. 

Both Ends Simply Supported 

Figure 23 plots the strains, e, , for the beam with both ends simply supported, calculated from 
equations (54) and (55) using the SPAR displacement outputs. The SPAR strain values presented 
in figure 23 were used to calculate the beam deflections, y., from deflection equation (12). 
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Table 12 presents these deflections and compares them to the deflections calculated from the SPAR 
program and those calculated from classical beam theory. 


Table 12. Comparison of deflections calculated from deflection equation (43) with those calculated 
from SPAR and classical beam theory; two-point supported aluminum tubular beam with both 
ends simply supported; P = -200 lb at beam center; n = 8. 





Deflection ( V] 

B ),in. 








B 

B 

B 

B 

B 

B 

B 


% 

3) 


3b 

3b 

(load point) 

3b 

3b 

3b 

3b 

SPAR 

0.0000 

-0.0333 

-0.0624 

-0.0834 

-0.0924 

-0.0834 

-0.0624 

-0.0333 

0.0000 

Classical 

beam theory 

0.0000 

-0.0320 

-0.0600 

-0.0798 

-0.0873 

-0.0798 

-0.0600 

-0.0320 

0.0000 

Deflection 

equation (43) 

0.0000 

-0.0318 

-0.0594 

-0.0787 

-0.0860 

-0.0787 

-0.0594 

-0.0318 

0.0000 

Difference 

from SPAR, 
percent 

Difference 

0.0000 

4.5045 

4.8077 

5.6355 

6.9264 

5.6355 

4.8077 

4.5045 

0.0000 

from classical 

beam theory, 
percent 

0.0000 

0.6250 

1.0000 

1.3784 

1.4891 

1.3784 

1.0000 

0.6250 

0.0000 


The deflection data listed in table 12 are plotted in figure 24 for visual comparison. Note that 
the deflection curve generated from equation (12) is very close to that generated from classical 
beam theory with a small difference of only 1.49 percent at the beam center. For the beam with 
both ends simply supported, which has no inflection points, the difference between the deflections 
calculated from equation (43) and those calculated from the SPAR program is only 6.93 percent at 
the beam center. Again, this deflection differential is caused by the shear effect, because when the 
shear modulus, G, was set to infinity (G — >°°), the SPAR program produced a deflection curve that 
was very close to that generated from classical beam theory. 

Plates 

This section presents the deflection predictions of plates under two different edge support 
conditions (four edges clamped and four edges simply supported). The deflections calculated from 
the Ko displacement theories are compared with those calculated from the SPAR finite-element 
analysis. 
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Four Edges Clamped 


Figure 25 shows the strains, e i , along the axis of symmetry of a clamped square plate (fig. 5), 
calculated from the SPAR displacement outputs using strain equations (54) and (55). The SPAR 
strain values presented in figure 25 were used to calculate the deflections, y t , along the axis of 
symmetry from cantilever beam deflection equation (12). Table 13 lists the results and compares 
them with the deflections calculated from the SPAR program. 


Table 13. Comparison of deflections along the axis of symmetry calculated from SPAR with those 
calculated from cantilever beam deflection equation (12); square plate with four edges clamped; 
P = -400 lb at plate center; n = 8. 


Deflection ( y t ), in. 



y 0 

A 

34 

34 

34 

(load point) 

34 

34 

34 

34 

SPAR 

0.0000 

-0.0067 

-0.0213 

-0.0380 

-0.0481 

-0.0380 

-0.0213 

-0.0067 

0.0000 

Deflection 
equation (12) 

0.0000 

-0.0066 

-0.0213 

-0.0378 

-0.0469 

-0.0367 

-0.0189 

-0.0031 

0.0047 

Difference, 

percent 

0.0000 

1.6442 

0.0939 

0.4474 

2.5572 

3.5263 

11.0798 

54.1106 



k 

Left-hand region 


>|<- — 


Right-hand region 

> 


Note from table 13 that the deflections calculated from equation (12) and those calculated 
from the SPAR program are very similar in the left-hand region (i = 0 - 4). Beyond the center 
loading point (i = 4), the differences between the deflections gradually increase in the right-hand 
region (i = 4 - 8), resulting in a small nonzero upward deflection, y 8 = 0.0047 in., at the right 
support. The reason for this nonzero upward deflection is the local strain perturbation created by 
the concentrated load at the station, i = 4 (plate center), which affects the deflection calculations 
for the right-hand region. For the present symmetrical loading case, the deflection data calculated 
from deflection equation (12) for the left-hand region can be used for the right-hand region to 
avoid error accumulation in the direct deflection calculations. 

When deflection equation (37) (corrected for the two-point supported beam) is used, the 
resulting deflection curve falls very close to the deflection curve calculated from the SPAR program 
for the right-hand region. As shown in table 14, however, the percent difference for the left-hand 
region increases slightly. 
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Table 14. Comparison of deflections along the axis of symmetry calculated from SPAR with those 
calculated from simple beam deflection equation (37); square plate with four edges clamped; 
P = -400 lb at plate center; n = 8. 


Deflection (»• ), in. 



y 0 

A 

3b 

3b 

34 

(load point) 

3b 

34 

3b 

3b 

SPAR 

0.0000 

-0.0067 

-0.0213 

-0.0381 

-0.0481 

-0.0380 

-0.0213 

-0.0067 

0.0000 

Deflection 
equation (37) 

0.0000 

-0.0072 

-0.0225 

-0.0396 

-0.0492 

-0.0396 

-0.0225 

-0.0072 

0.0000 

Difference, 

percent 

0.0000 

7.4627 

5.6338 

4.2105 

2.2869 

4.2105 

5.6338 

7.4627 

0.0000 


k 

Left-hand region 


>|< 


Right-hand region 

> 


The deflection data from tables 13 and 14 are plotted in figure 26 for graphical comparison. 
Note that the deflection curves generated from equation (12) and those generated from the SPAR 
program are very close in the left-hand region but gradually deviate as the right end is approached. 
The curve generated from equation (37) exhibits better overall predictions. 

Four Edges Simply Supported 

Figure 27 shows the strains, e, , along the axis of symmetry of a four-edge simply supported 
square plate (fig. 5), calculated from the SPAR displacement outputs using strain equations 
(54) and (55). The SPAR strain values presented in figure 27 were used as inputs to deflection 
equation (43) to calculate the deflections, y . , along the axis of symmetry. Table 15 compares these 
results to the deflections, y . , calculated from the SPAR program. 
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Table 15. Comparison of deflections along the axis of symmetry calculated from SPAR with 
those calculated from deflection equation (43); square plate with four edges simply supported; 
P = -400 lb at plate center; n = 8. 


Deflection ( yf ), in. 



yf 

yf 

yf 

yf 

B B 

34 34 

(load point) 

yf 

yf 

yf 

SPAR 

0.0000 

-0.0320 

-0.0622 

-0.0876 

-0.1005 -0.0876 

-0.0622 

-0.0320 

0.0000 

Deflection 

equation (43) 

0.0000 

-0.0330 

-0.0640 

-0.0898 

-0.1021 -0.0898 

-0.0640 

-0.0330 

0.0000 

Difference. 

percent 

0.0000 

3.1250 

2.8939 

2.5114 

1.5920 2.5114 

2.8939 

3.1250 

0.0000 


|< Left-hand region >l< Right-hand region > 


The deflection data from table 15 are plotted in figure 28 for visual comparison. Note that 
the deflection curve calculated from equation (43) and the deflection curve calculated from 
the SPAR program are quite close. The slight difference between the two curves is pictorially 
inconspicuous. 


DISCUSSION 

In-flight predictions of the deformed shape of the ultralightweight flying wing require the 
installation of multiple bending and distortion strain sensors, equally spaced on the surface of 
the wing tubular spar, for measuring the bending and distortion strains of the wings. The strain 
sensor data can then be input to theoretical displacement equations for calculations of slope, 
deflection, and twisting of the wing tubular spar at all strain- sensing stations, including the spar 
tip sensing station. 

Because the accuracy of the displacement (bending and torsion) equations have been validated 
by the classical beam theory and finite-element analysis, the equations are ready for application. 
To solidify confidence, however, the accuracy of the displacement equations must be checked 
extensively by testing several types of continuous and composite cantilever tubular beams under 
bending, torsion, and combined bending and torsion. The validation experiments are presently 
being conducted. 


When strain sensors are installed on the wing spar (or wing box), the displacement 
equations together with strain sensor data can be used to construct the in-flight deformed shape 
of an aircraft such as the Helios flying wing. Because they are lightweight, fiber optic sensors 
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may be the optimal choice for the embedded strain-sensing system on the wing spar surface. 

The use of a conventional strain gage sensing system is impractical, because it has numerous 
lead wires that add weight, thus it is too heavy for the lightweight flying wing to carry. 

CONCLUDING REMARKS 

Displacement equations were developed for a cantilever wing tubular spar under bending, 
torsion, and combined bending and torsion loading. These displacement equations were expressed 
in terms of strains measured at multiple sensing stations equally spaced on the surface of the wing 
spar. The principal results are as follows: 

1 . The displacement equations developed for the uniform cantilever beam were successfully 
validated for accuracy by classical beam theory and finite-element analysis. 

2. The deflections for the three types of cantilever tubular spars predicted from the 
deflection equations agreed quite well with those calculated from the finite-element 
analysis. 

3. The deflection equations developed for the uniform cantilever beam were extended to 
predict, with adequate accuracy, the deformed shapes of a two-point supported uniform 
beam (simple beam). 

4. By means of a multiple-line strain-sensing system, the deflection equations developed 
for simple beams can be applied to predict the deformed shapes of plates. 

5. Deflection equations were also developed for tapered, slightly tapered, and stepwise 
tapered cantilever beams. These displacement equations were applied to predict the 
deflections and twist angles of a tapered cantilever wing box (for which two strain- 
sensing lines are required) with reasonable accuracy compared to the finite-element 
predictions. 

6. The displacement equations and associated strain-sensing system (such as fiber optic 
sensors) create a powerful means for in-flight calculations of slopes, deflections, 
and cross-sectional twist angles of the wing spar at any strain-sensing station. These 
calculated displacements can be used to construct the deformed shape of the long-span 
flying wing. Ultimately, the calculated deformation data can be visually displayed for 
the ground-based pilot to monitor the in-flight deformed shape of the long-span flying 
wing. 
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Strain sensor 



X i - 1 X i 060098 


Figure 3. Instrumented cantilever tubular spar under combined bending and torsion. 



Figure 4. Wing box instrumented with a two-parallel-line strain sensor system. 




Figure 5. Flat panel instrumented with a multiple-parallel-line strain sensor system. 



/ =100.5 in. Nodes 3673 

a = 4 in. (Al, four-ply walls) Four-node elements 3636 

t = 0.02296 in. (Al, four-ply walls) Three-node elements 36 

a = 4.07749 in. (Helios spar caps) 

7 = 0.17784 in. (Helios spar caps) 

Figure 6. Helios tubular spar finite-element model subjected to a vertical tip load of P = 100 lb. 
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2 c n = 3 in. 

t = 0.02296 in. 

c i = c 0 ~ ( c 0 “ c n) [(*A/)//] 


Figure 7. Wing box finite-element model subjected to bending and torsion. 


Fixed 
clamped 
end 


Undeformed 


200 lb 



Free 

clamped 

end 


/ = 100.5 in. 
a = 4 in. 
t = 0.02296 in. 
c = a + 0.5 1 


Nodes 3673 

Four-node elements 3636 
Three-node elements 36 


060103 


(a) Both ends clamped. 


Figure 8. Finite-element model for an aluminum tubular beam subjected to a downward central 
vertical load of P = -200 lb. 


48 



Free simply 
supported 
end 


/ 


/ = 100.5 in. 
a = 4 in. 
t = 0.02296 in. 
c = a + 0.5 1 


Nodes 3673 

Four-node elements 3636 
Three-node elements 36 


060104 


(b) Left end clamped and right end simply supported. 



Free simply 
supported 
end 


/ 


/ = 100.5 in. 

Nodes 

3673 


a = 4 in. 

Four-node elements 

3636 


t = 0.02296 in. 

Three-node elements 

36 


c = a + 0.5 1 



060105 


(c) Both ends simply supported. 


Figure 8. Concluded. 
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w = 24 in. Four-node elements 2356 

h = 0.2 in. 

c = hi 2 = 0.1 in. 060106 

(a) Four edges clamped. 


I P = -400 lb at 



w = 24 in. Four-node elements 2356 

h = 0.2 in. Two-node elements 98 

c = hi 2 = 0.1 in. 000107 

(b) Four edges simply supported. 

Figure 9. Quarter panel finite-element model (with mirror image) for a square plate subjected to a 
downward central vertical load of P = -400 lb. 
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Figure 1 2. Bending strains, e ; , at different sensing stations, x t , calculated from the SPAR program; 
four-ply composite cantilever tube; n = 8; P = 100 lb. 



Strain-sensing station identification number, i 


(a) n = 2;P = 100 lb. 


Figure 13. Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (12); four-ply composite cantilever tube. 
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Deflection, V/, in. Deflection, v z -, in. 



Strain-sensing station identification number, i 


(b) n = 4;P = 100 lb. 



Strain-sensing station identification number, i 


(c) n = 8; P = 100 lb. 


Figure 13. Concluded. 
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Strain-sensing station identification number, i 


Figure 18. Comparison of twist angles, 0, , calculated from equation (36) using the SPAR 
program and deflections calculated from equations (35) (stepwise) and (30) (tapered); tapered 
cantilever wing box; n = 8; P = 100 lb; T = 50 x 17 in-lb. 


58 




012345678 


Strain-sensing station identification number, 


060123 


Figure 19. Bending strains, £ ( , at different sensing stations, x n calculated from the SPAR 
program; two-point supported beam with both ends clamped; n = 8; P = -200 lb. 



Strain-sensing station identification number, i 

060124 


Figure 20. Comparison of deflections calculated from the SPAR program, deflection 
equation (12), and classical beam theory; two-point supported beam with both ends clamped; 
n = S;P = -200 lb. 
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x 1CT 3 



Strain-sensing station identification number, i 

060125 

Figure 21. Bending strains, £ ( , at different sensing stations,^., calculated from the SPAR 
program; two-point supported beam with a clamped left end and simply supported right end; n = 8; 
P = -200 lb. 



Strain-sensing station identification number, i 

^ 060126 

Figure 22. Comparison of deflections calculated from the SPAR program, deflection 
equation (12), and classical beam theory; two-point supported beam with a clamped left end and 
simply supported right end; n = 8; P = -200 lb. 
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Strain-sensing station identification number, i 

060127 


Figure 23. Bending strains, £ ( , at different sensing stations, x ( , calculated from the SPAR 
program; two-point supported beam with both ends simply supported; n = 8; P = -200 lb. 
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SPAR 

Equation (12) 
Beam theory 


Figure 24. Comparison of deflections calculated from the SPAR program, deflection 
equation (43), and classical beam theory; two-point supported beam with both ends simply 
supported; n = 8; P = -200 lb. 
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x icr 3 



Strain-sensing station identification number, i 

060129 

Figure 25. Bending strains, e i , at different sensing stations, x i , calculated from the SPAR 
program; square plate with four edges clamped; n = 8; P = -400 lb. 



Strain-sensing station identification number, i 

J 060130 
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Figure 26. Comparison of deflections calculated from the SPAR program, equation (37), and 
deflection equation (12); square plate with four edges clamped; n = 8; P = -400 lb. 



x io~ 3 



Strain-sensing station identification number, i 

060131 

Figure 27. Bending strains, £, , at different strain-sensing stations, x t , calculated from the SPAR 
program; square plate with four edges simply supported; n = 8; P = -400 lb. 



Strain-sensing station identification number, i 

060132 


Figure 28. Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (43); square plate with four edges simply supported; n = 8; P = -400 lb. 
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APPENDIX A 
HELIOS FLYING WING 


The Helios prototype (fig. 1) is an ultralightweight, unmanned, solar-powered flying wing 
aircraft (1,600 lb) designed to fly at altitudes up to 100,000 ft. The Helios wing has 14 electric 
motors (2 hp each) that drive 79-in. -diameter composite propellers. The current necessary to power 
these motors and other systems is generated by high efficiency solar cells spread over the upper 
surface of the wing. The cruising speed of the Helios is in the range of 19 to 27 mi/h. 

The Helios wing has a chord length of 8 in. uniform over a 247-ft wingspan. The Helios 
wingspan is longer than the wingspan of the U.S. Air Force C-5 military transport (Lockheed 
Martin Aeronautics Company, Marietta, Georgia), which is 222 ft, and that of the Boeing 747 
commercial jetliner (Boeing Company, Seattle, Washington), which ranges from 195 to 215 ft 
depending on the model. The Helios wing has an area of 1,976 ft 2 and a maximum wing loading 
of 0.81 lb/ft 2 . 

The wing thickness is 11.5 in. from tip to tip, without taper or sweep. The outer wing panels 
have built-in 10-deg dihedral. The Helios is constructed of mostly composite materials such 
as carbon fiber, graphite epoxy, Kevlar® (E.I. duPont de Nemours and Company, Wilmington, 
Delaware), Styrofoam® (Dow Chemical Company, Midland, Michigan), and a thin, transparent 
plastic skin. The main tubular wing spar is fabricated with four plies (45/-45/-45/45) of carbon 
fiber composite. To increase bending stiffness, an additional 27 plies (0 g /90/0 8 /90/0 g ) of carbon 
fiber composite were added in the 60-deg regions of the top and bottom of the spar. The spar is 
wrapped with Kevlar and Nomex® (E.I. duPont de Nemours and Company, Wilmington, Delaware) 
for additional strength. The wing ribs are made of epoxy and carbon fiber. Shaped Styrofoam is 
used for the leading edge of the wing, and a durable clear plastic film covers the entire wing. 
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APPENDIX B 

INTEGRATIONS OF SLOPE AND DEFLECTION EQUATIONS FOR 

UNIFORM BEAMS 

This section explains the details of the integrations of slope equation (6) and deflection 
equation (9) for uniform beams. These details provide the steps that result in the final mathematical 
forms given by equations (7) and (10), respectively. 

Slope Equations 

The slope, tan 0(x ) , of the uniform beam in the region, x t , <x<x n between the two adjacent 
strain-sensing stations, [jq _ l ,x i } , can be obtained by integrating equation (3) with the constant of 
integration determined by enforcing the continuity of slope at the inboard adjacent strain-sensing 
station, x t , as 


tan0(x)= [ 

J x, 


d 2 y 
dx 1 


Slope increment 


dx+ tan0,._j 

Slope at X;_i 


£(X) 


dx + tan6 : 


i-1 


< X < jc f 


(B-l) 


The strain, e(x ) , in the region, x j _ l <x< x t , between the two adjacent strain- sensing stations, 
{*,_!,*, } , is assumed to be a linear function as 


£(*) = (£,■_! -£, ) 


x — X 

M 


7-1 


(B-2) 


Substitute equation (B-2) into equation (B-l), and carry out the integration as follows 
(ref. 4): 


1 r 

tan 0(x) = - 

n J . 


C yx i-i 

1 


-£,) : 


A / 


dx + tan 


£j_\(x Vj-i ) (£,_, £,) 


(x - v,_i ) 2 
2 A l 


(B-3) 


+ tan0 ( _. 


At the strain-sensing station, .y , one can write v, - .r, , = A / , and equation (B-3) yields the 
slope, tan 0\ = tan 0(x i ) ] , at the strain-sensing station,^. , as 


tan£f[= tan0(x,)] = - 

c 


£/_! (A/) -(£,._,-£,.) 


(A If 
2M 


+ tan 0 j x 


(B-4) 


After grouping the terms, equation (B-4) takes on the final form of the slope equation for the 
uniform beam as 
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tan0 ; =— (e ( ._, +e) + tan0 ; _, ; (i = 1, 2, 3, n) 

2c 


which is equation (7) in the text. 


(B-5) 


Deflection Equations 

The deflection, y(x), of the uniform beam in the region, * <x<x , between the two 
adjacent strain- sensing stations, , can be obtained by integrating slope equation (B-l) 

with the constant of integration determined by enforcing the continuity of deflection at the inboard 
adjacent strain-sensing station, x t _ { , as 

p x c x r* c x 

y(x) = tan 9(x)dx + y,_, = dxdx+ tan O j _ l dx+ y t _ x 

J Xj_[ 1 J X;_i J X,_j Q J Xj_i L — y 1 

v v ' Deflection v v ' v v J Deflection 

Integration of slope at x { _ x Deflection increment Deflection at x at Xj_ { 

due to tan0,_i 

Substitute equation (B-2) into equation (B-6), and carry out the integration as follows 
(ref. 4): 


T - i ^ x < x, 

(B-6) 


= 11.7 


e,_i -(£,_!-£, ) 


x - x 


i-i 


A / 


dx + tan 0 i , >dx + y ( , 


£ i-i( x x i -\ ) ( e ;- 1 e ;) 

(■ x-x ,_i) 2 


(x-x,._i ) 2 


■'i-i 




2 A l 

(x~Xj_ i ) 3 

6 A / 


+ tan# , \dx + y t _ 


(B-l) 


+ (x-x i _ l )tmd i _ l +y i _ l 


At the strain-sensing station, x t , one can writer, - x t _ x = A l , and equation (B-7) yields the 
deflection, y,[= >’(-*, ) I , at the strain-sensing station,^ , as 


= y(*i) = - 
c 


„ (A/) 2 N (A/) 3 

£, i (£, i -£,) 

2 11 ' 6 A/ 


+ (A/)tan6f_, + y ( _j 


(B-8) 


After grouping the terms, equation (B-8) takes on the final form of the deflection equation for 
the uniform beam as 


y> -i+ e ;) + A/ tan 0,._, +y,_j ; {i = 1, 2, 3, ...., n) 

6c 


which is equation (10) in the text. 


(B-9) 
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APPENDIX C 

DEFLECTION EQUATIONS WITH SHEAR EFFECTS 


For a long cantilever beam, the contribution of deflection caused by transverse-shear 
effect is usually very small. If the transverse-shear effect is desired, however, additional 
transverse- shear-induced deflection terms must be added to deflection equation (12). 

If Vj is the local shear force acting vertically at the strain-sensing station, x l , then the 
averaged shear strain, y i _ l , at the strain-sensing stations, x ( _, , induced by V i at the strain-sensing 
station, x t , is given by 


7 , 


AG 


(C-l) 


in which A is the beam cross-sectional area, and G is the shear modulus. 

The segment- wise transverse-shear-induced deflection at stations, relative to station x t , 
can be written as (A . The total transverse-shear-induced deflection, y s , at station x i is then 
obtained by summing up all the segment-wise deflections from the support station j = 0 up to 
j = (/' - 1 ). Namely, 


y'i = (A/)(7o + 7 1 + 72 + 7 3 + 74 + ■■■■ + 7,-1 ) (C- 2 ) 

The transverse-shear-induced deflection, y) , given in expression (C-2) can be added to the original 
cantilever beam deflection equation (12) to account for the shear effect: 


A/ 2 ^ \ 

>’i - G~\ 2e o + e \ ) 

6c 

_ A/ 2 , x 

y 2 - —— (5e 0 + 6e, + e 2 ) 
6c 


Al 2 r / , -i 

y 3 - — [8e 0 +6(2e! +e 2 ) + e 3 J 


y 4 [lle 0 + 6(3e, +2e 2 +e 3 ) + e 4 ] 


y. (transverse-shear effect term) 


+(A/)7 0 (C-3a) 

+(A/)(7 0 + 7j) (C-3b) 

+(A/)(7 0 + 7 1 + 7 2 ) (C-3c) 

+(A/)(7o + 7i +Y2+Y3 (C-3d) 
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y n = 


(A if 


6c 


n — 1 


(3 n - l)e 0 + 6^ (n - i)e i + e n 


+ (A /)(r 0 + 7, + 7 2 + 7 3 + 7 4 +.... + Y n -\ ) (C-3e) 


If the uniform cantilever beam is loaded by the vertical tip load, P, then the segment-wise 
shear strains will be constant along the beam span. Namely, 


/0=7l = /2 = /3 = /4= ••••=/„-! 


Then, in light of equation (C-4), deflection equation (C-3) becomes 


(C-4) 


A/ 2 . ■, 

>’i - ~Y~V- e o +e i) 

6c 

_ A/2 (e A \ 

y 2 - —— (5e 0 + 6e, + e 2 ) 
6c 


A/ r / , -i 

y 3 - — [8e 0 +6(2e, +e 2 ) + e 3 J 

y 4 =-^— [lle 0 + 6(3e, +2e 2 +e 3 ) + e 4 ] 


y. (transverse-shear effect term) 


+(A/)7o 


+2(A/)7 0 

+3(A/)7 0 

+4(A/)7 0 


(C-5a) 

(C-5b) 

(C-5c) 

(C-5d) 


(A/) 2 

6c 


(3 n - l)e 0 + 6 - i)e i + e n 


+n(Al)y 0 


(C-5e) 
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APPENDIX D 

INTEGRATIONS OF SLOPE AND DEFLECTION EQUATIONS FOR 

TAPERED BEAMS 

This section explains the details of the integrations of slope equation (22) and deflection 
equation (24) for tapered beams. These details provide the steps that result in the final mathematical 
forms given by equations (23) and (25), respectively. 

Slope Equations 

The slope, tan 0{x ) , of the tapered beam in the region, x ( ._, < x < x,. , between the two adjacent 
strain-sensing stations, {x^, ,x ( } , can be obtained by integrating equation (3) (c is replaced with 
c(x) ) with the constant of integration determined by enforcing the continuity of slope at the inboard 
adjacent strain-sensing station, x t _ x , as 

rx d 2 y cx e(x ) 

tan 0(x) = — T dx+ tan6>. , = t/x + tan# , ; x,,<x<x, 

Jx i->dx~ v— ^ j, hc(x) (D-l) 

v -v ' Slope at X{_i 

Slope increment 

In the region, x j _ l < x < x ( , between the two adjacent strain-sensing stations, {x, p x ( } , the 
strain, e(x) , and beam half depth, c(x) , are assumed to be of the following linear functions: 

£(x) = - (£,_, - £,) ' V ^ 1 (D-2) 


c(x) = c i _ l -(c i _ l -c i )^-^ (D-3) 

Substitute equations (D-2) and (D-3) into equation (D-l), and carry out the integration as 
follows (ref. 4): 
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tan(x) = [ 

J X: 


, X- X. 1 

A/ 


C/-i - (c,--i - C i) 


dx + tan 6 


x - Xf 

M 


c i-i( e /-i £/) e ,-i( c i-i C () 


(c f _i — c,.) 


(x-x ; _,) + 


A/ 


A/ 


(A/) 2 


log 


M 


(a--x ( _,) + c,_ 


-logc M l + tan0, 


i-1 


e '-' £/ - (X - X ■ , ) + A/ e 7 |C; ^ <! log 




£,r 


( c ,--i ~ c ,-) 

M 


(x — X ( ._, ) + C ; _j 


-logc M I+tan0 f 


i-1 


(D-4) 


At the strain-sensing station, x ( , one can write x ; - x ( . , = A/ , and equation (D-4) yields the 
slope, tan 0 ( .[= tan 0(x j ) ] , at the strain-sensing station, x ( . , as 


tan 0 i - tan 0(x l ) = 


£ j _ l -£ i . . . . £;_,C ; - e.c 


A/ + A/ ^ (log c,. - log C ( _j ) + tan 0._, (D-5) 

(Cj-l-C,-) (C/-1-C/) 


After grouping the term, one obtains the final form of the slope equation for tapered 
beams as 


tan 0 i = A l 




£/-l C ! 


£,C ; _l , C; 


(^-1 “ C i f 


-log — 


+ tan0 


( i = 1, 2, 3, n) (D-6) 


which is equation (23) in the text. 


Deflection Equations 

The deflection, y(x) , of the tapered beam in the region, x ; , <x<x,, between the two 
adjacent strain-sensing stations, {x ; _,,x,}, can be obtained by integrating slope equation (D-l) 
with the constant of integration determined by enforcing the continuity of deflection at the inboard 
adjacent strain-sensing station, x ( _j , as 


f x C x C x r x 

y(x) = tan 0(x)dx+ y i _ l = — —dxdx+ tan6 j _ 1 dx + 

j X: , L — 1 j X: 1 j Xt 1 n( J X: 1 


Integration of slope 


Deflection 

at X:_, 


c(c) 


y t - 


; A. 1 - ■* - x, 


Deflection increment 


Deflection at x 
due to tan0,_! 


Deflection 

at X:_ 1 


(D-7) 
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Substitute equations (D-2) and (D-3) into equation (D-7), and carry out the integration as 
follows (ref. 4): 


X— X; 


rx rx £ i 1 ( £ '~ 1 £ ') A / rx 

y(x) = I I v dx+ I tan d^dx + y^ 

J X ;_ , J X:_t X X , * X:_, 


=r 

•'-'Vi 


C: | -(C, | -C;) — 

.-1 A/ 

* I £,_! -e, . . . , - e,c, ! 

' 1 ' (jc — ) + A/ 1 1 ' ' 1 1 


( c ,--i “ c /) 




log I '— ^(* - Vi) + C t _ x 1 - logc,.., 


A/ 


+ (* _ V_i ) tan 6 j , + y ( _. 


e /-l C ; £ / C (-l 


£ ‘ 1 £ ' -(X- JC;_j ) 2 + A l 


(Ci-l-C,) 


2 r 




( c ,--i -Cf) 

M 


~i C i- i-C/) 


A/ 


U- V- I ) + C i - 


x- 


lo §l (6 '~ a/ Cf ) (* ~ V-i ) + c,-- 1 j - lQ g V-i 


(^■-i 

M 




+(x-x i _ 1 )tmO i _ l +y i _ 1 


(D-8) 


At the strain-sensing station, x , one can write x - x ; _, = A/ , and equation (D-8) yields the 
deflection y.[= y(x)] at the strain-sensing station, x t , as 


y, = y(x/) = g ' (A/) 2 - (A/) 2 £/ - |C '' £, ' Cf - 1 


2 (^-i - <■'« ) 

- c /) 


(Cl-l-Ci) 


X 


A/ + C- | 


M 

+ A/tan0 ; _ 1 +y ; ._. 


lo§ (C 'a/ C ) a/ + C; - 1 ) " log c, -‘ 


+ (£ £z 1 Z5) A / 

M 


= (A/) 2 <! £ '"‘ _ £ ' - £ ' _ £ ^3~’ [c, (logc, - logc w ) + (c,._j - c,.)]| + A/ tan0 ( _, 

[ZtC ( | C,) (C ; _J c t ) J 

which can be written in a final form of the deflection equation for tapered beams as 
V = (A/) 


+ V-i 


(D-9) 


£ i-l - £ i 

£ i-l C i £ i C i-\ r 

2(c,_i - c ( ) 

7' 

1 

jri 

U) 

1 


c.-iog — + (c,-_i-c i ) 
c. , 


+ A/ tan#,. t + y ( ._j 
(z = 1, 2, 3, n) 


(D-10) 


which is equation (25) in the text. 
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APPENDIX E 

EXPANSIONS OF LOGARITHMIC TERMS 


When the cantilever beam is slightly tapered (that is, c ( — > c, , ), the logarithmic terms in 
slope equation (23) and deflection equation (25) will approach zero [that is, log(c,/c / _ 1 ) -> 0]. 
Therefore, the perturbation method must be used to expand the logarithmic term, log(c, /c,_ 1 ) , in 
the neighborhood of (c ; /c,_,) ~ 1 to yield nonzero mathematical expressions, which are derived in 
this section. 


Expansion of log(c ; /c ;1 ) for Slope Equation 

When the beam is slightly tapered (that is, c, — > c ; _, ), the logarithmic term log (c j /c j _ l ) can be 
expanded in series in the neighborhood of c ( ./c,. , ~ 1 (that is, log(c ; /c ; _ 1 ) ~ 0 ) (ref. 4) as follows: 


log— - 
C,_, 







c ) ~ E- I , i ) 2 

2 c i-i 3c?_j 



>1 (E-l) 


Retaining the terms up to the second order terms in (c^, - c,) 2 , the logarithmic term log(c, /c i _ 1 ) 
becomes the following form when c, approaches c) , (c, — > c [ _ ] ): 



c 


;-i 


c i - 1 


c, 


2c, 2 




c, 


>1 


c 


i-1 


(E-2) 


Equation (E-2) is identical to equation (26) in the text. 

Expansion of [trlog^/c^) + (c M - c, )] for Deflection Equation 

When the beam is slightly tapered (that is, c, — > c, , ), the term, [c, log(c ( /c, ,) + (c, , - c ; )] , 
in deflection equation (25) can be expanded in the neighborhood of (c,/c,_,)~l (that is, 
log(c ; /c,_,) ~ 0 ). In light of the logic, /c, ,) expansion equation (E-l), the expansion of 
[c, log(c,/c, ,) + (c, , - c,)] can be carried out up to the third order terms in (c, , - c,) 3 when c, 
approaches c,_, (c,. — > c, , ) as follows (ref. 4): 
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(E-3) 


Neglecting the forth order term in - c, ) 4 , one obtains the following final expression for 
[c, logfc, /6' _| ) + (c,_| -c,)] when c t approaches c, , (c, -> c, , ) : 


c. 


log— + (c,_, -c,) 


c . 


/-I 


(Cj-l -Cjf 

6c ,-l 


3 + 


(Cj-1-C,-) 


"i-1 


— -> 1 (E-4) 

c ,--i 


which is equation (29) in the text. 
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